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The “Palace of Discovery” 
at the Paris Exposition of 1937 


By JESSE W. M. DUMOND 


California Institute of Technology, Pasadena, California 


URING the summer of 1937 it was my 

pleasure and privilege to visit on many 
occasions the ‘‘Palace of Discovery”’ of the Paris 
Exposition. This vast, splendid and instructive 
exhibit of scientific principles and progress was 
essentially the idea of that ‘grand old man”’ of 
French science, Jean Perrin. The details were 
admirably worked out with the wholehearted 
cooperation of many of the most important 
French scientific workers who gave freely of 
their time and energy. 

The purposes of M. Perrin in organizing this 
scientific exhibition are so well expressed in his 
own words that I merely quote him 
here in my own translation. 

“It should be made evident to 
all that the recent and prodigious 
progress of our civilization has had 
its source in pure and disinterested 
research directed toward the dis- 
covery of the unknown. Thus, no 
one should leave this exhibit with- 
out realizing that no part, I re- 
peat absolutely no part, of the 
vast electrical industry of today 





were it not for the discoveries of 
the electrical currents and their properties by 
the Italian, Volta; the Frenchman, Ampere; and 
the Englishman, Faraday. 


would exist 


“We wish at the same time to disseminate the 
idea that this country, that any country, has a 
primordial and very practical interest, however 
great its financial difficulties may seem, in giving 
to its first class research workers the sums of 
money, relatively insignificant as they are com- 
pared to a national budget, which they need for 
their research instead of economizing absurdly 
on the seed with which future harvests are sown. 

“Finally pinning our hopes on 
the plain people, who in every 
country have shown in the most 
touching and moving manner their 
faith in the value of Science, we 
wish to stress among the adoles- 
cents who come with eager curi- 
ousity to see our experiments, the 
opportunities of vocational choice 
which may lead to the repetition 
of such miracles as the rise of 
Faraday from the position of book- 
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Fic. 1. Mounting of cathode ray to show sound waves. 


binder-workman to that of the most eminent 
physicist of his country. 

“The ‘Palace of Discovery’ built for this 
triple purpose was organized by French scien- 
tists. Those who visit it will see, however, that 
it is as international as is the science which it 
summarizes. Save for any forgetfulness or 
omissions which it shall be our constant effort 
to correct, we have tried to display the work of 
the great research men of all nations.”’ 

M. Perrin goes on from here to point out the 
solidarity of effort in world research which 
becomes even more apparent when we consider 
present living workers. For example, a discovery 
in nuclear chemistry may frequently be the 
result of the conscious collaboration of researchers 
in far distant lands all attracted at the same 
time toward the same enigma. He concludes by 
expressing the hope that a crusade in favor of 
scientific research such as the recent one in 
Belgium may spread to all countries in order 
that they may reap both the ideal intellectual 
and natural benefits that are sure to follow and 
he closes with the sentence, ‘““‘Do you not see 
that then Peace will follow of necessity ?”’ 

The Science Exhibit was located in the “‘Grand 
Palais,’’ the huge building which has in the past 
always been used for the annual exhibitions of 
French painting and sculpture. It was divided 
into sections of mathematics, physics, chemistry, 
astronomy and the biological sciences. It is 
impossible in a short article to do justice to all 
of these different sections. The excellent French 
journal La Nature in its issues of the first and 
fifteenth of September, 1937, has given fairly 
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complete and accurate descriptions accompanied 
with many photographs and diagrams of the 
exhibits and experiments in the sections of 
physics and chemistry, respectively. I wish here 
to speak especially of the exhibits classified under 
physics. Taking for example the sections of 
mechanics and structure of matter, these were 
divided into sub-sections with the following 
names: Mechanics or Hall of Galileo (M. A. 
Marcellin), Crystallography (Mauguin and 
Wyart), States of Matter (J. Basset), Molecular 
Magnitudes (N. Choucroun), Fluorescence and 
Phosphorescence (F. Perrin, Rouault and 
Coustal), Negative Electrons and X-Rays (H. 
Hulubei, G. Y. Cauchois), Cosmic Rays (P. 
Auger and P. Ehrenfest), Radioactivity (M. De- 
bierne), Atomic Synthesis (F. Joliot). I have 
purposely given after the title of each section 
the name of the French physicist who devoted 
time and energy to make it a success although 
these names were nowhere in evidence in the 
exhibit. It seems to me that it was this whole- 
hearted and self-sacrificing devotion and coop- 
eration of the French representatives of pure 
science and research which is the most instructive 
object lesson of the exhibit. Their work would of 
course have been sadly handicapped had they 
not been generously supported with funds and 
technical assistance. 

One could, of course, go on to mention all the 
other divisions of physics, such as geometrical 
optics (Fabry, Bayle, Arnulf, Guerian, Pencio- 
lelli), physical optics and the velocity of light 
(Cotton, Andant), the Hall of Ampere and 
Faraday (Cotton, Tsai, Leprince-Ringuet, Poir- 
son); oscillatory phenomena (Langevin, Lucas), 





go photographs in this paper 

and the one used on the cover 
were supplied by the editors of the 
French periodical *‘La Nature,”’ 
and we take this opportunity to 
express our appreciation for their 
kindness and cooperation.—THE 
EDITOR 
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physics of the globe (Maurin, Salles) with their 
yarious subdivisions but space does not permit. 
The place was a veritable physicists’ paradise 

a show appealing at once to the most highly 
sophisticated savant and to the general public 
there was a wealth for everybody. I shall mention 
therefore just a few exhibits and experiments 
taken at random. 





Fic. 2. General view of the ‘Hall of Negative Electrons’”’ 
with the Lenard ray tube in the immediate foreground 
operated by the 250 kv 30 milliampere constant tension 
generator seen above. Fast electrons issue from the thin 
window in the bottom end of the tube L where they pro- 
voke luminescence in the air and in minerals which are 
placed below and which exhibit superb colors. Booths 
containing automatically operated exhibits can be seen 
lining either wall of this large room. 


In the Hall of Galileo not far from the Foucault 
pendulum is a circular table 
turning night and day. Small plumb bobs and 
vessels of fluid placed at different radial distances 


large rotating 


from its axis of rotation showed that a citizen of 
this little rotating world would at each position 
form a quite different impression of the vertical 
and horizontal from that formed by the sta- 
tionary observers in the room. The lesson was 
driven home with curious vitality by the fact 
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Fic. 3. Section of the large display of crystal models. 


that grass and other plants planted in earth dis- 
tributed all over the surface of the rotating disk 
had grown in the direction of the plumb bobs and 
normal to the water surfaces but inclined relative 
to the earth’s vertical by greater angles the 
farther the plants and grass were from the center. 
The effort of growth of the plant was always 
exerted against what for it was the prevailing 
gravitational field. Thus was both centrifugal 
force and Einstein's principle of equivalence 
driven home. 

I cannot forget the demonstrations of change 
of state, the critical point, the continuous passage 
from gaseous to liquid states, etc. so ably ex- 
plained to the spectators. Nor can I pass over 
the elegant microprojections permitting anyone 
to see the Brownian movement, discontinuous 
thickness of thin films, etc., and to understand 
through clear explanations the evidence for the 
reality of the molecular structure of matter. The 
entrancing new fields of monomolecular layers 
dimensional fluids furnished some 
delightful experiments. Such layers apparently 


or two 


can behave as two-dimensional solids (as Devaux 
has shown) or as two-dimensional liquids or gases 
(Marcellin, Langmuir, Blodgett). 

Under crystallography a bewildering multitude 
of models made to a uniform scale of 400 million 
times magnification showed how x-rays reveal 
teach us 


the structure of materials and most 


valuable practical information as to why, for 
example, diamond is hard and graphite is soft 
though they are made of exactly the same atoms. 


The value of this discovery in pure science for 
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the understanding of nature, even in biology or 
in the practical field of the invention of new 
materials with desired physical properties is a 











Fic. 4. Wilson cloud expansion chamber to demonstrate 
ranges of alpha-rays. Built by Surrugue. 


lesson whose value we think should sink deeply 
into the minds of those who hold the purse 
strings of potential funds for research. 

The highly practical sections on electrolysis 
and electrophoresis as well as the section on 
photoluminescence can only be mentioned with 
warm praise for lack of space. However, the 
Hall of Negative Electrons was of such exciting 
interest to me, I must tell about it in some 
detail. Mlle. Cauchois and M. Hulubei are to be 
congratulated for the splendid results of their 
really very ambitious efforts. Nearly all the 
exhibits which lined the walls in booths on either 
side of a large room were automatic so that the 
visitor, by pushing appropriate buttons and 
levers and reading the accompanying printed 
explanations, could instruct himself as to the 
nature and laws of thermionic emission, the 
deflection of electron beams by electric and 
magnetic fields, the ways of determining the 
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specific charge em, photoelectric emission, 
electric discharges in different gases at different 
pressures, X-ray generation, x-ray absorption, 
photographic effects of x-rays, ionization of gases 
by x-rays, and many other experiments jn 
generous profusion. Especial mention should be 
made also of the contributions to this exhibit 
from the work of Manne Siegbahn and of the 
electron diffraction tube operating automati- 
cally which J. J. Trillat of Besancon set 
up. The latter was arranged so that by push- 
ing a button anyone could see as often as he 
liked a beam of electrons diffracted by a crystal 
forming beautiful clear and regular ‘‘Laue’”’ spots 
on a fluorescent screen. The important results 
which Cauchois and Hulubei have achieved in 
the development of the curved crystal focusing 
x-ray spectrograph which has recently permitted 
the discovery of element 87 by these workers 
were exhibited along with many examples and 
forms of this powerful new tool. As a climax 
there was in the center a large high voltage 
cathode-ray tube fitted with a Lenard window 
from which could be made to issue into the air 
fast electrons which falling on an assortment of 
minerals made the latter emit a gorgeous display 
of colors. 

For those who read without running printed 
explanations and spectra displayed about the 
“Hall of Negative Electrons’’ made evident such 
highlights of the history of atomic physics during 





Fic. 5. Trillat’s permanently evacuated and automati- 
cally operated electron diffraction tube. The fluorescent 
screen at the left end of the tube showed the von Laue 
patterns of the electron waves diffracted by a crystalline 
lattice very clearly to all spectators who pressed the 
button. 


JOURNAL OF APPLIED_PHYSICS 


























the twentieth century as the Bohr atom model, 
the concept of atomic number of Moseley, the 
wave particle duality of electrons and of radi- 
ation based on the immortal experimental find- 
ings of Gray, Compton, Davisson and Germer, 
and G. P. Thomson together with its enunci- 
ation in Heisenberg’s principle of indetermina- 
tion. 

One is torn between the limitations of space 
and the aversion for omitting to speak of even 
one of so many well conceived displays; the 
section on radioactivity and atomic synthesis, 
for example, with its Wilson cloud chambers in 
operation and with the original apparatus of the 
Curies; the section on cosmic rays with the 
automatically started cloud chamber set off by 
counters so skillfully realized by Auger and 
Ehrenfest; the projection views accompanied 
with automatic phonograph transcriptions ex- 
plaining the epoch making discoveries and 
photographs taken with the Wilson chamber by 





Fic. 6. High luminosity spectrographs and some models of 
organic chemicals from the chemistry section, 


workers in this field from America, Engiand, 
France and many other countries with a fine, 
generous impartiality. There was the huge elec- 
trostatic belt generator of Van de Graaff type 
which formed the ‘‘piéce de resistance,”’ as one 
entered the great central rotunda of the exhibit 
which so impressed the French public that it 
Was very amusingly used in political cartoons 
in the newspapers in which the debates in the 
Luxembourg over financial policy were repre- 
sented as the crashing disruptive discharges of 
this machine. 
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Fic. 7. The constituents of air including the rare gases 
were separated and exhibited by means of electric dis- 
charge tubes continually for every spectator to see in this 
apparatus. Section on chemistry. 


By the time one had reached the section on 
optics one was prepared for anything and again 
one was not disappointed. There were beautiful 
and lucidly clear explanations of geometrical 
optics showing the rays going through lenses or 
reflected by mirrors or showing the formation of 
the mirage or of the rainbow; the experiment of 
the measuring of the velocity of light with mirrors 
turning on a jet of air and a beam only thirty 
meters long; every type of interference and dif- 
fraction experiment of which it seemed to me I 
had ever heard and many others all in perfect 
adjustment and clearly explained so as to make 
the wave properties of light self-evident, and a 
host of other beautiful experiments making clear 
the scattering of light, the colors of the sunset, 
the polarization and rotatory polarization of 
light, double refraction and its multitudinous 
applications, and so on almost it seemed without 
end. 

The Ampere-Faraday hall contained a huge so- 
called ‘“‘Shomopolar’’ generator which I believe 
was said to furnish 50,000 amperes in some 
experiments. These tremendous currents flow- 
ing through huge conductors were made to 
perform weird tricks at various places in the 
room and visitors were warned to leave their 
watches at a distance. Measurement of the mag- 
netic fields with the ‘‘Cotton balance,’ a clear 
demonstration of the Zeeman effect, experiments 
on the magnetic birefringence of liquids, on 
falling bodies in a magnetic field and a very large 
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Fic. 8. Cauchois curved crystal spectrograph for demon- 
strating absorption and emission x-ray spectra. G, x-ray 
tube; O, optical bench; P, photographic plate holder. 
From the ‘Hall of Negative Electrons.”’ 


magnetic Wilson cloud chamber of de Broglie and 
Leprince-Ringuet were a few of the other ex- 
hibits which used the source of the immense 
electrical current just mentioned. 

Of the sections on sound and ultrasonics and 
other vibratory phenomena we have not room 
to speak nor can we say but a word about the 
there were 
the 
cellular formation of thermoconvective vortices. 


fine section on meteorology where 


striking experimental demonstrations of 
And we must not forget that my description 
the than 


the great display for instance of the 


says nothing of subdivisions other 
physics; 
progress in biology, medicine, surgery, heredity, 
microbiology, etc., ete., nor can I speak of the 
subdivisions of astronomy nor of mathematics 
with its museum of precise mathematical instru- 
ments or its little round treasure room whose 
frieze was decorated with the value of + com- 
puted to heaven only knows how many decimals 
spiralling again and again around. And last but 
least 


not the subdivisions on chemistry and 


industrial chemistry so vital to our modern 
civilization must be passed undescribed. 

But I have said enough perhaps to point out 
the moral to this tale which many a reader may 
already have guessed. In this work-a-day world 
I think recently one can detect a misunder- 
standing which may have grown out of the too 
speedy and careless popularization of science—a 


misunderstanding of the nature, value, and role, 
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as Perrin put it, ‘‘of pure and disinterested research 
directed toward the discovery of the unknown.” 
There is no doubting the meaning of his well 
chosen words. As our own Jewett has said ‘‘Yoy 
can’t have applied research until you have built 
up by pure research a body of knowledge to 
apply.”’ Experience has shown (and an exhibit 
such as the French ‘Palace of Discovery”’ is a 
this) that these 
frontiers of the known which must be ceaselessly 
pushed forward in order that all our technique 
of civilization may flourish and grow in its wake 
can only be pushed forward by disinterested 
research. The immediate practical so-called utili- 
tarian goal must therefore not be our exclusive 
aim or we risk nipping in the bud the Newtons 
and the Faradays and the Maxwells and the 
Diracs and the Heisenbergs. It is the lesson of 


marvelous way of teaching 


their part in the world picture that such a 
science exhibition should aim to teach. Whether 
we believe that government (as in France) or 
private initiative (as in America) is the best 
agency to encourage pure research is certainly 
a question which I believe is still open for debate 
but it is certain to me in either case that the 
value of ‘‘pure and disinterested research directed 
toward the discovery of the unknown” must be 
driven home to the public unless civilization is 
to die a slow death from malnutrition in these 
United States. The physicists of America have 
two opportunities to do their part in this cause 
one is the exhibition in San 
Francisco, the other the exhibition in New York. 


in the near future 


Let us hope they will see the need and rise to 
the occasion. 








Fic. 9. Reflection, refraction, and grating diffractions ol 
short radio waves. 
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Reviews of Recent Research 


Study of Gamma-Radiation from Boron 
Bombarded by Protons 


HE knowledge of the properties of the ex- 

cited states of atoms and molecules secured 
by experiment forms a large part of the basis of 
the present theory of atomic and molecular struc- 
ture. It is now believed that experimental data 
on the properties of nuclear states of excitation 
will aid in developing and testing an adequate 
picture of nuclear structure. One method of de- 
termining the properties of excited nuclear states 
is by studying the intensity and energy of the 
gamma-radiation emitted in transitions to lower 
lying states. Those transmutations involving the 
capture of fast protons by light nuclei are espe- 
cially useful in this connection as they result in 
the production of excited states of only one 
nucleus or at most two if the bombarded element 
has two isotopes. In various laboratories, studies 
of the gamma-radiation transmutations 
have been attempted by securing cloud chamber 


from 


photographs of secondaries (Compton recoil elec- 





An electron-positron pair produced in the gas of a cloud 
chamber by a gamma-ray of 16.1 Mev energy produced 
in turn by bombarding boron with fast protons. The 
electron member of the pair bending toward the inter- 
section at the two stereoscopic views of the cloud chamber 
has 9.5 Mev kinetic energy while the positron has 5.6 Mev 
kinetic energy. The production of the pair required 1.0 
Mev of the gamma-ray's energy. The curvature is due to 
a magnetic field of 2500 gauss. 
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trons and electron-positron pairs) produced in 
thin sheets of material such as lead, carbon, or 
aluminum by the radiation in question. 

Recently* W. A. Fowler, E. R. Gaerttner, and 
C. C. Lauritsen of the W. K. Kellogg Radiation 
Laboratory of the California Institute of Tech- 
nology have found in this way that when boron is 
bombarded with protons, gamma-rays of energy 
16.6+0.6 Mev, 11.8+0.5 Mev, and 4.3+0.3 Mev 
with relative intensity 1/7 :1:1 are emitted. 
The nucleus formed in this case is an excited 
state of carbon which can decay to normal car- 
bon directly or in a double transition through a 
state at 4.3 Mev above normal carbon. This state 
has been found to be produced in numerous other 
transmutations. 


* 


Penetrating Power of High Velocity 
Electrons Studied 


LECTRONS fast travel 130 
meters in air can now be produced in 
abundance by means of artificial nuclear disinte- 


enough to 


gration processes. These serve to fill in the gap 
between the extreme energies found in cosmic 
rays (which easily traverse the entire atmosphere 
of the earth) and the much lower energies of 
the naturally occurring radioactive substances. 
Experimenters can therefore make measurements 
of the phenomena associated with the penetra- 
tion of electrons through matter over the entire 
spectrum of energies. 

During the past two years detailed measure- 
ments have been carried out at the University 
of Michigan to determine the ability of electrons 
of energies from 0.5 to 17 million electron volts 
to penetrate sheets of lead and carbon. This 
energy range is especially interesting because it 
includes a transition from one type of absorption 
to another. Theory predicts the following: At 
the lowest energies the loss of energy suffered 


* Phys. Rev. 53, 628 (1938). 
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by an electron is due almost entirely to the 
production of ions along its path or to collisions 
with other electrons. At very high energies the 
latter is entirely overshadowed by losses due to 
the production of x-rays or high energy quanta. 
When lead is used as the absorbing material 
the kinds of 
absorption falls within the energy range studied: 


transition between these two 
for example, at 3 Mev the loss is due 75 percent 
Mev it 
30 Mev it is due 75 percent to radiation. The 
the the Michigan 
experiments, the latest of which are reported by 
A. J. Ruhlig and H. R. Crane in the April 15th 


to ionization, at 10 


is 50-50, and at 


values for losses found in 


issue of the Physical Review, are higher than the 
theoretical predictions over the entire range of 
energy studied. They run from 3 to 4 times the 
theoretical loss at 0.5 


\lev to about 1.5 times 


the theoretical loss at 15 Mev, and from the 
trend of the values it appears that as the energy 
is increased the agreement between experiment 


and theory will become better. One of the 





Cloud chamber photograph showing the change in 
curvature of an electron path before and after passing 
through a sheet of lead. Electrons enter the chamber from 


the lower edge of the picture. 
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difficulties in the interpretation of the exper. 
that the pursue 
crooked paths inside the absorbing material, 
making it 


mental results is electrons 
hard to estimate the true distance 
traveled. This is particularly true for low cnergy 
electrons, and may very well account for at 
least some of the apparent excess in the energy 
loss. 


The 


figure will give an idea as to the method of 


cloud chamber picture shown in the 
experiment. Electrons, whose paths are curved 
field, entering the 
chamber from the lower edge of the picture, 


by a magnetic are seen 
Some of them pass through the sheet of lead, 


which appears as a horizontal line across the 
center. The curvatures of the paths before and 
after passing through the lead tell how much 
energy had with and 
how much it lost in traversing the lead. Many 
times backward from the 


lead. This is due to one or more large deflections 


each electron to begin 


electrons “‘bounce”’ 
caused by close approaches of the electron to 
lead nuclei. The two electrons entering 2nd and 
3rd the left in the 
about 180 degrees with very small loss in energy. 


from fivure are deflected 


The electrons shown in the figure have energies 


of the order of 4 Mev. 
* 


Proton Induced Radioactivities 
HILE 


operation in this country have confined 
the 
produced by high speed deuterons and alpha- 


most of the cyclotrons now in 


their attention to nuclear disintegrations 


particles, the one at the University of Rochester* 
has been devoted to the investigation of reactions 


produced by high energy protons. Using a 


proton beam whose energy is 4 Mev (‘which has 
more recently been increased to 5.5 Mev) a 
large number of radioactive isotopes have been 
produced. The particular feature of this work is 
that a previously unknown type of reaction has 
been discovered in which a proton entering a 
nucleus ejects from it a neutron. In almost all 
the resulting nucleus is radioactive, in 


Cases 


some cases emitting positrons and in other 


cases, electrons. In this paper results are reported 
on the radioactivities observed in proton bom- 


* Described by L. A. 


DuBridge, ef a]. Phys. Rev. 53, 
447-454 (1938). 
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The cvelotron of the University of Rochester with a 
27-inch chamber. It is rated to produce 7 Mev protons. 


bardment of the following elements: N, O, Cr, 
Co, Zn, Se, Mo, Pd, Cd, In. Particularly strong 
radioactivities were observed in targets of oxygen 
and selenium and the energy relations involved 
in these reactions have been studied in some 
detail. 

The positron emitting isotopes formed by the 
proton-neutron reaction are of particular interest. 
In this case after decay the radioactive isotope 
returns to the initial isotope. The mass relations 
are very simple and lead to the prediction that 
the minimum proton energy at which the reaction 
may take place is equal to 1.8 Mev plus the 
maximum the emitted 
This relation has been experimentally verified. 

* 
Electrostatic Generator Operating Under 
High Air Pressure 
A RECENT paper in the Physical Review 
describes improvements in the high voltage 


generator developed at the University of Wis- 


energy of positrons. 


consin for experiments in nuclear physics. The 


the and 


generator is of electrostatic type, 
' PD. B. Parkinson, R. G. Herb, E. J. Bernet and J. L. 


McKibben, Phys. Rev. 53, 642 (1938). 
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operates in a horizontal cylindrical tank 20 feet 


long and 55 feet in diameter, under an air 


pressure of 100 Ib. in, this pressure giving a 
maximum voltage four times as high as that 
obtained with air at atmospheric pressure. The 
addition of a few kilograms of CCl, or CClhF» 
vapor to the air in the tank has been found to 
increase the maximum usable voltage from 2.1 
to 2.4 million volts; so that the upper limit to 
the voltage is now set by a breakdown along 
the accelerating tube, the charging belts, or the 
radial 
sparking across the 18-inch air gap between the 
high potential electrode end the wall of the tank. 
(senerator 


supporting insulators rather than’ by 


voltage can be easily set at any 
desired value up to the maximum and is held at 
a given value to within less than 1 percent by 
an automatic voltage stabilizer. 

Two balloon fabric charging belts running 
along the axis of the tank from one end feed 
charge to the high potential electrode which is 
located near the center of the tank. Extending 
from this electrode along the axis of the tank to 
the opposite end is the evacuated accelerating 
tube, which consists of a number of corrugated 
porcelain cylinders separated by metal elec- 
trodes. Protons from the high potential electrode 
travel down the accelerating tube, and pass out 
through the evacuated tube at the end of the 
tank, where they are at present being used for 
proton-proton scattering experiments, and for 
measurements on nuclear transmutations caused 
by proton bombardment. 





The Wisconsin electrostatic generator. 





Recent Progress in Supersonics : 


By WILLIAM T. RICHARDS 


Rockefeller Institute for Medical Research 
New York, New York 


MUST begin with an apology. In going over 

the development of supersonics during the 
past five years I found that it is quite impossible, 
in a short address, to give credit where credit is 
due in all cases. What I am about to say, there- 
fore, will be an attempt to trace the development 
of several main lines of investigation, and much 
interesting work which is not intimately asso- 
ciated with these main lines cannot be men- 
tioned. I hope that any of you who consider that 
I have sinned either by omission or by commis- 
sion will make your views known to me as loudly 
and as emphatically as possible. In an expansion 
of this address, which will later appear in the 
Reviews of Modern Physics, 1 shall give a more 
scholarly exposition. I shall speak first of the 
propagation of sound in gases, liquids and solids. 
Later I shall devote a few paragraphs to the 
chemical and physical effects of sound of high 
amplitude. 

The theory of the propagation of sound in dis- 
persive gases has been very fully treated by 
numerous writers. In general they follow either 
the statistical or the thermodynamic discipline. 
In general also they consider the dispersion due 
to the failure of a part of the heat capacity to 
follow the adiabatic cycle of the sound wave. To 
illustrate this as specifically as possible let us 
consider chlorine, a gas which is partly excited 
vibrationally at room temperature and which 
has, in consequence, a heat capacity at constant 
volume greater than 5 2 of R per mole. If the 
velocity of sound is measured at low frequencies, 
the rate at which the small pressure impulses of 
the sound waves succeed each other is less than 
the rate at which translational energy is con- 

*An Address to the Acoustical Society of America at 


Ann Arbor, November 30, 1937. Read by the Society's 
president for the author. 
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verted into vibrational energy, and the effective 
heat capacity is the equilibrium or ‘‘zero fre- 
quency” heat capacity. In the ‘‘zero frequency” 
range the gas is dispersive only because of second 
order effects such as viscosity, heat conduction, or 
finite amplitude of the waves. When the fre- 
quency has been raised to within one order of 
magnitude of the rate of conversion of transla- 
tional into vibrational energy, dispersion from 
this cause becomes manifest. The reason for the 
dispersion is the failure of some of the molecules 
of the gas to alter their vibrational energy to con- 
form to the rapidly changing local temperature. 
This results in a diminution of the effective heat 
capacity. The pressure-temperature cycle of the 
sound wave remains adiabatic, but is no longer at 
maximum entropy. A phase angle has appeared 
between the propagation of the pressure and the 
density waves, and the tangent of this angle gives, 
of course, the absorption coefficient. At circular 
frequencies one order of magnitude or more 
above the rate of conversion of translational into 
vibrational energy, no appreciable effect of the 
vibrational energy remains. The velocity of sound 
has reached a new and higher value, which re- 
mains constant even at very high frequencies, 
and corresponds to a heat capacity of 5 2 R for 
chlorine. Both the thermodynamic and the sta- 
tistical treatment agree in describing this disper- 
sion by the methods of modern chemical reaction 
kinetics, although this is not necessarily inherent 
in the thermodynamic viewpoint. Hence they ar- 
rive at essentially similar conclusions, which are 
currently termed ‘‘the collision theory of acous- 
tical dispersion.” 

Bourgin is responsible almost single-handed for 
the development of the statistical theory. In 
1928 he wrote the only complete theory of dis- 
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persion for waves of infinitesimal amplitude that 
I have seen.' He considered viscosity, quantiza- 
tion of the collision process, radiation from ex- 
cited molecules, and dispersive gaseous mixtures. 
He followed this by expanding his treatment of 
gaseous mixtures, and, in 1932, wrote the first 
theory of the three-state gas.” This is a gas having 
translational energy and two sets of internal en- 
ergy states which are characterized by different 
rates of adjustment. In chlorine the two states 
would be the vibrational and the rotational en- 
ergy. Finally, he has recently compared the 
statistical and the thermodynamic methods.’ 
Considering the brilliance of this achievement 
Bourgin should be honored first among those 
writing on dispersion, yet not one experimenter 
has, so far as I am aware, taken advantage of this 
work. I think that this is only because the sta- 
tistical method involves reasoning of considerably 
greater mathematical difficulty than the thermo- 
dynamic. We experimenters are simple folk, 
easily awed by a row of summation signs, and we 
have chosen the simplest language we could find. 
Professor Bourgin has had to wait ten years for 
the recognition which is his due. He may have to 
wait another ten. But since a statistical method is 
the most logical way of describing a statistical 
process, he must ultimately triumph. 

The thermodynamic method was first pub- 
lished in articulate form by Einstein in a paper of 
great beauty and simplicity dealing with disper- 
sion in a dissociating gas. Kneser, ten years 
later, applied the same method to a pressure- 
density phase lag in a nondissociating gas.® 
Kneser made a serious error which Rutgers, 
leaning even more heavily on the exposition of 
Einstein, corrected. I found it necessary to ex- 
tend Einstein’s treatment to a three-state gas in 
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Richardson's hot-wire probe apparatus for measuring 
acoustic absorption in gases. 
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order to prove that it was the dissociation reac- 
tion and not the heat capacity which caused dis- 
persion in nitrogen tetroxide.’ Luck extended 
Einstein's treatment to describe absorption in a 
real dissociating gas,* and I followed his lead by 
correcting for van der Waals forces in nondisso- 
ciating gases and gaseous mixtures.’ I also worked 
out the behavior of four- and five-state gases, but 
this, judging by the profound apathy with which 
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Knudsen’s arrangement for measuring the absorption of 
sound in a gas by the intensity method. 


the publication’? was greeted, was carrying 
matters too far. The work of Landau and Teller'' 
makes it appear that such monstrosities will not 
be necessary. The method of Einstein—and as a 
reviewer I must insist that it be called the method 
of Einstein and not the Kneser-Rutgers method 

has yet to be corrected for viscésity and heat 
conduction. 

The very stimulating contribution of Herzfeld 
and Rice,” which was prior to all the work men- 
tioned above save only that of Einstein, does not 
fall readily into either of these two classifica- 
tions, although it also postulates a_ collision 
mechanism. 

It should be apparent from what I have said 
that a believer in the collision theory may find 
appropriate expressions to describe his measure- 
ments in either Bourgin’s or Einstein's language. 
In both sufficiently explicit expressions have been 
written to enable him to identify the cause of any 
dispersion or absorption he may find, provided 
only that his experimental technique is equal to 
the occasion. 

Fortunately, perhaps, for the vitality of the 
subject, some investigators, who have obtained 
results which they believe inconsistent with the 
collision theory, proceed to deny its validity. I 
think especially of Richardson" and of Penman" 
in England. I should be entirely in sympathy 
with this enterprise if I had not personally en- 
countered far more startling irregularities which 
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[he diffraction of light by sound waves. The upper 
figure shows the intensity distribution with progressive 
waves and the lower with standing waves. (After R. Bar. 


later proved to be apparatus effects. I shall men- 
tion one or two of these in a few minutes. A col- 
lective work entitled Aberrations of Sonic Inter- 
ferometers would have a useful place in acoustical 
literature. In view of this I shall here take the 
high-handed stand that measurements which do 
not agree with the collision theory are experi- 
mentally unsound. If this stand can be refuted, 
my interest in the subject will be greatly re- 
freshed. 

Before leaving the theoretical aspect of the 
propagation of sound in gases I should like to call 
attention to a tendency on the part of experi- 
menters to overidealize the situation. For ex- 
ample, expressions developed for ideal gases are 
often cheerfully used over wide temperature and 
pressure ranges. The most common of these 
abuses is the so-called “‘reduction to zero degrees 
centigrade” by treating the velocity as if it were a 
linear function of the square root of the absolute 
temperature. For carbon dioxide this is not even 
a close approximation. Indeed, ‘‘zero frequency” 
velocity measurements appear to offer an accu- 
rate means of evaluating the second virial coeftf- 
cient. Again, there is the tendency to regard the 
waves generated by a piezoelectric or magneto- 
strictive oscillator as plane, purely sinosoidal, and 
of infinitesimal amplitude. They certainly have 
none of these characteristics to any high degree of 
approximation. Fay, in a valuable and little-read 
paper,” has shown that, for waves of even very 
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moderate amplitude, the attenuation due to the 
conditions of energy transfer may be twenty 
times that due to viscous loss. This attenuation is 
due solely to the fact that the sound wave is a 
sound wave, and that the condensed portions are 
traveling faster than the rarefied portions in pro.- 
portionality to the amount of matter in each. 
Pielemeier has concerned himself with the experi- 
mental aspects of this problem.'® Further, we 
tend to neglect the effect of viscosity and heat 
conduction. This effect is small except in tubes 
having a diameter of the order of magnitude of 
the wave-length, but we should not forget that 
it is present. In cases like oxygen at room tem- 
perature, where the low and high frequency 
velocities differ by only a few tenths of a percent, 
it may be necessary to reckon with it. Finally, the 
possibility of radiation from a dipole as a means 
of energy loss must not be overlooked. In carbon 
dioxide the radiation must be slow, and has been 
shown experimentally not to affect the acoustical 
measurements. In other gases, however, which 
have large permanent dipoles and dispersive re- 
gions in low frequency ranges, it could well 
appear. 

There is little to be said concerning the theory 
of sound propagation in liquids and solids. In 
liquids, the frequency of molecular collisions, if 
one can use such a term in a semi-rigid lattice 
structure, is so great that there is little chance of 
detecting dispersion due to heat capacity in fre- 
quency regions at present available for experi- 
ment. Moreover, if dispersion were found, its 
interpretation would present grave difficulties 
since the kinetic theory of liquids, though defi- 
nitely conceived by the efforts of Fowler, Eyring 
and Bernal, is still in the foetal stage and cannot 
be expected to give assistance to other branches 
of the science. It appears that, conversely, 
measurements on the scattering of sound in 
liquids will lead to a substantiation of some of 
these theories, as Lucas!’ has recently suggested. 
Dissociating liquids, such as acetic acid, might be 
expected to show dispersion owing to the chem- 
ical reaction involved. Also there should be dis- 
persion, closely analogous to the anomalous dis- 
persion of radio waves, during the liquid-glass 
“transition.”’ Similarly, in some crystals, where 
free rotation of certain groups begins below the 
melting-point, measurements of the dispersion of 
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sound should provide interesting information 
concerning the rate of energy transfer to these 
groups. I am not aware, however, that any of 
these phenomena have been reported. 

By far the most popular instrument for meas- 
uring the velocity and absorption of sound is the 
sonic interferometer introduced by Pierce.'> It 
has been modified for liquids and solids, and has 
been used to frequencies of half a megacycle 
with magnostrictive and twenty megacycles with 
piezoelectric oscillators. Nearly a hundred papers 
have appeared dealing with various aspects of this 
method alone. As evidence has accumulated it 
has become increasingly clear that it is not suffi- 
cient to set up an oscillator, calibrate it, and 
measure the position of points of maximum 
electrical reaction to secure numbers representing 
the velocity and absorption of sound which have 
absolute numerical significance. The wave pat- 
tern given off by the oscillator is not a pure plane 
sinosoid, since various parts of the surface are 
not vibrating in phase with each other. Grossmann 
has shown that, in consequence of this, a plane 
reflector will give maximum reaction on the 
oscillator at a position which is not exactly an 
even number of half-wave-lengths from the sur- 
face of the oscillator, but is a compromise with 
the various components of the wave pattern.!® 
In addition, all parts of the oscillator are sending 
out waves, and, if measurements are made in a 
confined space, wall reflections occur. When all 
this is considered it is clear that very complicated 
reaction-distance patterns may be found, and 














_Diffraction of a supersonic wave by a wire grating. The 
piezoelectric source is on the left. The first-order diffraction 
spectrum is in a southeast direction, while the second-order 
spectrum is in a southwest direction. (After R. Bair. 
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that they must be interpreted only with the 
greatest caution. Observers are often deceived by 
a reproducibility of the order of 0.001 percent 
into a belief that the numerical accuracy is 
equally great. Another difficulty which is too 














Diffraction and reflection of a supersonic wave in the 
surface between p-xylol (at the top) and water (at the 
bottom). Rays are incident from the northwest and are 
reflected to the northeast. Waves parallel to and just below 
the interface are going to the right. (After R. Bir.) 


little appreciated is the effect of varying pressure 
on the frequency of an oscillator. It does not ap- 
pear reasonable that a change from 250 mm to 
500 mm pressure of acetic acid vapor on a mas- 
sive nickel tube a foot long should change its 
frequency by 300 cycles, yet I have observed this 
phenomenon repeatedly and reproducibly. Cases 
are on record where no greater increases of pres- 
sure than this have caused a high frequency crys- 
tal wholly to stop oscillating. No doubt the effect 
is a multiple one, involving not only the direct 
damping of the oscillator, but also adsorption 
phenomena of a chemical nature which alter the 
relation of the crystal and its mount. These diffi- 
culties notwithstanding, the sonic interferometer 
is an instrument of high precision. It is capable of 
giving absolute velocities to better than 0.1 per- 
cent, and absorption coefficients, owing especially 
to the examination of the reflection coefficients 
involved by Hubbard,’ to about 10 percent. 
There are three welcome newcomers among 
the instruments suitable for measuring the veloc- 
ity and absorption of sound in various media. 
One of these, which was devised by Debye and 
Sears,”! employs the diffraction of light by the 
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Visibility of 
The frequency 
Hiedemann and Ch 


stationary ultrasonic waves in m-xvlene. 
is 5 megacyvcles. This picture is by E. 
Bachem 


alternate condensation and rarefaction of matter 
in the sound wave. Since standing waves are not 
required, the method is very simple, and it has 
been widely employed, especially by investigators 
interested in liquids. It is steadily undergoing 
improvement from a technical standpoint. The 
theory of the diffraction of light by sound waves 
has, in addition, attracted wide attention among 
mathematical physicists, the interpretations of 
Brillouin” Nath”™ being 
favored by most writers. The accuracy of this 


and of Raman and 
method appears to be almost as great as that 
obtainable with the Hubbard and Loomis” inter- 
ferometer. It has the additional advantage of 
requiring small volumes of material; Bar has 
measured the velocity with only 1 cc of D.O.” 

A second new method is that using a movable 
hot wire probe which indicates regions of maxi- 
mum and minimum motion by the temperature 
of the wire. This was devised by Richardson" and 
his co-workers, who recommend it especially for 
absorption measurements in gases. A thermopyle 
has also been used in this way. A systematic com- 
parison of the performance of this instrument 
with the interferometer would be valuable. Since 
it has not been made, and since the hot wire 
probe has been used by few investigators, it is 
difficult to estimate its accuracy. A paper which 
has recently appeared in this journal describing 
an investigation of the nodal points in wind 
instruments by the hot wire method shows that 
it may be valuable in several branches of 
acoustics. 

The most recent method depends upon the 
direct visibility of ultrasonic waves of even very 
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moderate amplitude. With intensities as low as 
2 watts per square centimeter standing waves of 
sound are visible to the eye in transformer ojl. 
With the ‘‘Schlieren-Methode” they may be 
photographed at even smaller intensities jn 
liquids, glasses, transparent crystals, and even in 
gases.”° Such photographs may be used to deter- 
mine the velocity of sound in conjunction with a 
microphotometer and a knowledge of the dimen- 
sions of the apparatus, gained either by direct 
measurement, or by comparison with a similar 
photograph of a substance of known properties, 





The formation of a fog at liquid-air interface by intense 
high frequency sound waves. (Photograph by Wood and 
Loomis. ) 


The velocity of sound may in this way be meas- 
ured in a hundredth of a second or less, which 
should be invaluable for a substance like formic 
acid vapor, which undergoes rapid thermal de- 
composition at temperatures where its vapor 
pressure is great enough to permit acoustical 
measurements. Bergmann, Heidemann, and 
others have recently developed interesting stro- 
boscopic variations of this procedure. The im- 
portant feature of these methods is that they 
sidestep any difficulties due to irregularities of 
the wave form. 

The 


brilliant success by Knudsen?’ for determining 


reverberation method, used with such 
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absorption coethcients at audible frequencies has 
been extended into the ultra-audible range.** It 
would be aesthetically pleasing if a method could 
be found to measure the group velocity, since 
this would give the dispersion directly. 

Turning now to the results of measurements, 
we must first consider the information to be ob- 
tained from the so-called ‘“‘zero frequency” or 
equilibrium velocities. Gases and vapors, if a 
suitable compensation is made for the van der 
Waals forces, yield in this region data which 
should be a valuable confirmation of the analysis 
of band spectra. For carbon dioxide and carbon 
disulfide such confirmation is trivial, but for 
gases like ethylene, where the analysis is difficult, 
accurate measurements of the velocity of sound 
at low frequencies should be really helpful in 
view of the difficulty of determining the heat 
capacity by other methods. The spectral analysis 
dioxide, chlorine, 
oxygen, and several other gases, has been checked 


of carbon dioxide, sulfur 
in this way. In liquids, since the heat capacity 
cannot at present be evaluated theoretically, the 
significance of the velocity of sound is necessarily 
more superficial, and in fact appears to be valu- 
able only for the determination of the adiabatic 
compressibility coefficient in nonelectrolytes. In 
electrolytes, since accelerated charges are in- 


volved, information concerning the masses of 


ions and their valences may also be obtained. 
Parthasarathy has been very industrious in 
measurements on a large number of organic non- 
electrolytes in the hope of establishing relations 
between the velocity of sound and molecular 
constitution, but has not had conspicuous suc- 
cess." The data on solids are meager. 

By far the most attention has been focused on 
the dispersion and absorption of sound in poly- 
atomic gases, since these data throw new and 
valuable light on the transfer of energy by im- 
pact. Too much has been published on this sub- 
ject to permit its detailed discussion in the time 
at my disposal. There are, however, certain con- 
clusions which result from the work of Knudsen, 
Kneser, Eucken, Pielemeier, Hubbard and their 
associates, and of Mr. Reid in my laboratory. | 
shall have to state these conclusions somewhat 
too categorically, but it is my belief that they are 
essentially correct. 
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(1) Dispersion and absorption of sound are to 
be expected in pressure, temperature, and fre- 
quency ranges at present accessible to experi- 
ment in any polyatomic gas in which an appre- 
ciable amount of vibrational energy is excited. 
The reason for this is primarily the inevitable 
inefficiency of the conversion of translational 
energy into vibrational energy by collision. In the 
vast majority of cases the effect of a collision will 
be merely to knock the molecules apart without 
bringing about any conversion at all. For the 
collision must convert one quantum of energy 
into vibration, not merely to perturb one of the 
two molecules sufficiently to set its atoms vibrat- 
ing very slightly. Furthermore, a steric factor 
must enter, since a diatomic molecule must be 
hit on its end to produce vibration. These are 
crude statements, but I am told that a more de- 
tailed study of the collision process produces a 
very similar picture. 

(2) The effect of increasing pressure is in ac- 
with the the kinetic 
theory of heat. An increase of pressure causes 


cordance indications of 
more collisions per unit time, and consequently 
pushes the dispersive region to higher acoustical 
frequencies. 

(3) The effect of increasing temperature is to 
increase exponentially the rate of conversion of 
translational energy into vibrational energy. The 
proportionality is not, however, that given by 
the Boltzmann factor for the particular energy 
state involved, but indicates that a larger energy 
than hy is required. The energy in excess of hy is 





Pattern produced by the reflection of ultrasonic rays by 
concave mirror. After E. Hiedemann and Ch. Bachem. 
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roughly analogous to the activation energy of a 
chemical reaction, but it has been shown thoret- 
ically to depend on temperature, so the analogy 
cannot be taken seriously. Landau and Teller 
show that only molecules of high relative trans- 
lational energy can effect conversion. Due to the 
phenomenon of ‘trace catalysis,’ which I am 
about to discuss, these temperature coefficients 
are very difficult to determine with accuracy. For 
this reason they are different, at present, in 
Germany from what they are in this country. I 
do not think that this should cause either sur- 
prise or distress. 

(4) Owing to the pioneer investigations of 
Knudsen, which have been amply confirmed by 
others, it appears that very small traces of foreign 
gases may have enormous effect on the efficiency 
of the collision process, and consequently on the 
behavior of a gas in the dispersive region. The 
effect of the addition of these traces is always to 
shift the dispersive region to higher frequencies. 
Any other result would necessitate a revision of 
the collision theory. Kneser and Knudsen*® in- 
terpret this result, I think correctly, as meaning 
that a collision of a molecule of the dispersive 
gas with a molecules of the “trace catalyst” is 
more effective in exciting the vibration of the 
dispersive gas than is a collision with one of its 
own species. The effect of changing concentration 
is not uniform, sometimes indicating that bi- 
molecular collisions are involved, and sometimes 
trimolecular, and so on. Eucken*! prefers to call 
this chemical excitation as opposed to physical 
excitation by impact with nonreactive gases such 
as He. Now phenomena like this have been 
familiar to chemists for a long time, and have 
distressed us greatly. If we are to think about it 
at all we must obviously go deeper than con- 
siderations of mere impact, since the average 
kinetic energy of all the molecules of the gas is 
the same. It is fairly reasonable to turn to rota- 
tional energy, and examine the molecules with 
small moments of inertia and consequently large 
rotational quanta. Chief among these is hydro- 
gen, which appears to be an effective catalyst, 
and deuterium, which is only about one-eighth as 
effective as hydrogen with ethylene. Water and 
ammonia, and perhaps even methane, appear 
promising candidates, and do not disappoint us 
in experiments, although their interrelationships 
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are not as simple as we could wish. But when 
Knudsen finds that ethyl alcohol is more effective 
with oxygen than any of a large number of trace 
catalysts he has tried, it becomes necessary to 
admit that rotational energy cannot explain, 
everything. For if, as Kistiakowsky and others 
have recently shown, there is no free rotation 
about the carbon —carbon bond in ethane, it js 
scarcely conceivable that there should be free 
rotation in ethanol, and the moment of inertia of 
the molecule, must be, therefore, very great. It 
seems to me, then, that we must turn in despera- 
tion to some idea of excitation by resonance be- 
tween the sum of a set of frequencies in ethanol 
and the vibrational frequency of oxygen, for it is 
only in this way that I can envisage the adding 
together of several small quanta to make a large 
one. The torsional quanta of ethanol must be 
considerable, and these may participate in the 
exchange. During the collision the energy in a 
number of degrees of freedom in ethanol must 
flow over to the oxygen atom. A statement like 
this makes me uncomfortable, but its equivalent 
has been in use for many years in chemical reac- 
tion kinetics. 

(5) Noinformation has been obtained concern- 
ing the rates of adjustment of the rotational 
energy or of the van der Waals forces, since both 
appear to follow the acoustical cycle to the high- 
est frequencies yet employed. 

(6) Chemical .reaction rates may be success- 
fully determined by the acoustical method, but 
it seems that only gases with low heat capacities 
are suitable for investigation at present. Indeed 
the dissociation of nitrogen tetroxide is the only 
success to date, and this is a gas which is excep- 
tionally easy to purify and manipulate, and hasa 
low heat capacity as well. 

I shall close with a brief survey of the effects of 
intense sonic and supersonic radiation. This 
branch of the subject, which originated with the 
work of Wood and Loomis,” is in a somewhat 
insecure position because no compelling reason 
has been found for many of the phenomena which 
occur. Nevertheless, a host of publications have 
appeared describing these effects, and attempting 
to give them some theoretical explanation. It 
soon became apparent that emulsification, the 
destruction of bacteria, the freeing of antibodies 
from bacterial cells, the acceleration of chemical 
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reactions, the production of luminescence and 
several other phenomena were dependent for 
their existence upon cavitation, which is the 
production of bubbles of gas in a liquid.** In gas- 
free liquids, or in liquids under high hydrostatic 
pressures, most of these disappear, and others, 
like emulsification, proceed only feebly and under 
restricted conditions. I have been told by every 
mathematical physicist I know that the analysis 
of cavitation is a task beyond the ability of pres- 
ent day mathematics. It is therefore necessary to 
examine the role of cavitation in the light of 
what experiment shows. 

We may profitably fix our attention on the 
acceleration of chemical reactions because this 
usually requires a definite quantity of energy, 
which any proposed mechanism involving cavita- 
tion must be able to supply. The most common 
type of reaction to be reported is of the oxidation- 
reduction type, for which the production of 
jodine from sodium iodide may serve as an ex- 
ample. Since it is possible to measure the in- 
tensity of sound quantitatively,* it is possible to 
show by straightforward thermodynamical rea- 
soning that no sound waves have ever been pro- 
duced with sufficient pressure amplitude to affect 





Jets of spray being sent out from a tube in which the 
walls are vibrating as stationary waves. Such vibration was 
produced by Wood and Loomis with high frequency sound 
Waves of great intensity. . 
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The propagation of ultrasonic energy through a dia- 
phragm having a diameter equal to 8/10 wave-length. 
Picture on the left shows the waves while the one on the 
right shows the acoustic rays. (After E. Hiedemann and 
K. Osterhammel. 


a chemical reaction appreciably owing to the 
temperature difference between the compressed 
and rarefied portions of the waves. It is possible, 
perhaps, that the high temperatures developed 
at the surface of cavitated bubbles, which Mari- 
nesco has shown to reach 250°C,** may locally 
accelerate reactions to an extent which appears 
as a measurable acceleration of the reaction in 
the liquid as a whole. This might account for the 
acceleration of reactions which are not of the 
oxidation-reduction type, such as hydrolyses. 
But it has repeatedly been observed that the 
irradiation with intense sound of air-saturated 
water produces hydrogen peroxide if cavitation 
is permitted, but has no result if it is suppressed. 
Here we are not dealing with the hastening of a 
process which leads to equilibrium, but with the 
creation of a compound which requires the ex- 
penditure of some 3 electron volts, and no tem- 
perature rise due to bubble-collapse, bubble- 
heating and the like can bring it about by any 
mechanism known to me. It appears, however, 
that cavitation is accompanied by severe elec- 
trical disturbances. Just as the production of a 
spray of liquid by an 


‘ 


‘atomizer’ causes a separa- 
tion of charges, the production of a spray of holes, 
which is cavitation, has the same result. Mr. 
Hobart and I have shown, in a paper to be sub- 
mitted to the Journal of the American Chemical 
Society, that large random electrical disturbances 
may be detected in a cavitating liquid. Over the 
distances comparable to molecular dimensions it 
appears reasonable to believe that these may be 
great enough to produce ionization, combination, 
and so on. Our proof is unfortunately not com- 
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plete, but it is clear that a ballo-electrical ex- 
planation of the role of cavitation in accelerating 
reactions is at present tenable. It also supplies a 
the 
liquids which Chambers has studied. 


welcome explanation of luminescence of 

There are many industrial applications of in- 
tense sound waves ranging from the destruction of 
fogs to the ageing of whiskey. Manufacturers of 
cosmetics and of photographic plates find intense 
sound useful in preparing stable and homogene- 
ous emulsions. The milk people are interested 
both from the standpoint of homogenizing and 
of bacterial destruction. Paints and varnishes 
have had a crack at it. In fact, about 1932 there 


was a feeling in the air that anyone who manu- 





factured anything, with the possible exception of 
horn buttons, was either installing a supersonic 
outfit or wishing that he had the money for one. 
The chief beneficiaries of this movement were the 


electric power companies. In some cases, it js 
true, better emulsions are produced by sound 
waves than by any other known method. The 
coagulation of fogs by airplanes zooming about 
with supersonic generators may have a great 
future. The smelting industry may find sound- 
waves indispensable for the degassing of metal 
melts. But the electrical production of sound 
waves is appallingly wasteful, and for the most 
part I believe that they will be supplanted by 
more efficient mechanical devices. 
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The Study of Extreme Pressures 
and their Importance in the 
Investigation of Engineering Problems 


BY THOMAS C. POULTER 


Executive Director, Research Foundation of Armour Institute of Technology, Chicago, Illinois 


ERHAPS the earliest 
high pressure work 
of note is the cele- 
brated attempt of the 
Florentine Academy to 
whether 


find water is 


compressible. Using a 
lead container in which 
they probably had much 
than a_ thousand 
pounds they 
concluded that water was 
an incompressible liquid. In 1762-64 Candon 
carried out some experiments in which he proved 


less 


pressure, 


that water was compressible, and in spite of the 
fact that water is sometimes referred to today 
as an incompressible liquid, experiments have 
been carried out in which water has been com- 
pressed to one-half its normal value. The fifty- 
year period ending 1869 was one in which a great 
deal was accomplished in the field of high 
pressure in spite of the fact that very few good 
numerical results were obtained. Since then a 
great deal of activity has been shown in the 
field of extreme pressure investigation, and the 
literature contains hundreds of articles published 
by a great many investigators so that at the 
present time there is a wealth of high pressure 
data available. 

If we examine the curve which represents the 
advance that has been made in the maximum 
pressures employed in this work in the past two 
hundred years, Fig. 1, we find that it is rising 
at such an increasing rate that it makes even the 
high pressure investigator wonder what the 
future may lead to in this field. In spite of the 
fact that man has until recently been unable to 
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develop pressures as high as one million pounds 
per square inch, such pressures are frequently 
occurring about us in our everyday experiences. 
They are, of course, confined to an extremely 
small space, and have a very short duration. 
Pressures of a few hundred thousand pounds per 
square inch are occasionally produced in the 
hypoid gear in the normal operation of an 
automobile, in roller or ball bearings, under the 
wheel of a glass cutter, the bullet of a high 
powered rifle striking a solid object, and in many 
other cases where two hard objects strike each 
other. Such extreme liquid or gas pressures are 
of course much less common, but even gas pres- 
sures of nearly a million pounds per square inch 
are probably occurring every day as particles of 
meteoric matter traveling with a velocity of 
thirty miles per second, collide with the very 
light upper layer of the earth’s atmosphere. 
Recently an experiment was carried out in 
the high pressure laboratory of the Research 
Foundation of Armour Institute of Technology 
in which a pressure of one million five hundred 
thousand pounds per square inch was developed. 
This was the highest pressure ever produced 
experimentally in high pressure investigations. 
The pressure at the center of the earth is esti- 
mated to be about two hundred 
thousand atmospheres, or only thirty-two times 


three million 


as great as the maximum pressure that has been 
produced in the laboratory. It is, of course, 
experimentally impossible to produce in the 
laboratory pressures corresponding to those at 
the center of the earth, but this increase in range 
of pressures may clear up many puzzling ques- 
tions which have formerly been considered im- 
possible of solution. Very little attempt has been 
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made in the past to utilize extreme pressure 
experiments in the solution of everyday prob- 
lems, but they are rapidly assuming a place of 
increasing importance, and in such experiments 
lie the answers to many questions of greatest 
commercial importance. 
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A brief survey of the work that has been 
studied by means of the high pressure equipment 
now in the 
Institute of 


the high pressure laboratory of 
Armour 
Technology will serve to show the great diver- 


Research Foundation of 
sified application of such investigations. The 
work was started in an attempt to study a 
reported effect of pressure upon the chemical 
reaction of sulphuric acid on zine to form zine 
sulphate and liberate hydrogen.' This investiga- 
tion was carried to a pressure of thirty thousand 
atmospheres, during which it was found that the 
effect we set out to study did not exist. In the 
course of the investigation we found many inter- 
esting things including a little known chemical 
reaction of hydrogen on sulphuric acid to give 
hydrogen sulfide and water. Other chemical reac- 
tions were investigated under extreme pressures, 
such as the hydrolysis of sucrose. The rate of 
this reaction was found to decrease with increase 
in pressure. From the information gained as a 
result of our previous work, new pressure equip- 
ment was developed in which it was possible to 
carry on experiments in cylinders with glass, 
quartz or diamond windows,?~* which will with- 
stand pressures of one-half million pounds per 
square inch. In connection with this development 
many interesting phenomena were observed, 
such as the ability of glass or quartz windows 
one-quarter inch thick to be bent to a radius of 
curvature of four inches as many as ten times 
before being broken. The development of such 
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windows made it possible to study the effect of 
pressures upon the optical rotation of optically 
active compounds, which effect was to increase 
the rotation of all compounds studied, some of 
them being increased to threefold their norma] 
Both and levorotatory com- 
pounds were investigated. The windows also 


value. dextro- 
made it possible to study the phosphorescence of 
zine sulfide,» which effect is to decrease the 
intensity of the phosphorescence by a factor of 
one-half for a pressure of thirty thousand atmos- 
pheres. It was also found that a rapid change in 
pressure would produce a bright glow of the 
zinc sulfide. Aside from the intensity, the fluo- 
rescent and phosphorescent properties of zinc 
sulfide were very little affected by extremely 
high pressures. It, therefore, provided a con- 
venient method of studying the effect of extreme 
pressures upon the radioactive decomposition$ 
of various materials by actually observing scintil- 
lation. These measurements confirmed the nega- 
tive results of similar investigations by other 
observers using quite different methods. Further 
investigations made possible by the pressure 
windows were the determination of the effect of 
pressure upon the index of refraction® of certain 
liquids. This that the 
Lorenz-Lorentz relation those 


investigation showed 
holds for sub- 
stances investigated. The quartz windows were 
further utilized in making absorption spectra 
measurements of the effect of pressure upon the 
compressibility of the neodymium atoms at the 
various The effect of 
pressure upon living organisms has been studied, 
with the that thousand 
atmospheres are necessary to kill bacteria, but 


electron energy levels. 


observation twelve 
as the complexity of the living organism increases 
the pressure necessary to destroy life decreases. 
Hydra and planaria were found to withstand 
pressures of from ten to twenty thousand pounds 
per square inch without any serious damage. The 
effect of high pressures was to precipitate some 
of the colloidal constituents of the organism, and 
a study was made of the precipitation of other 
sulfur, gold, 
hydroxide, molybdenum blue, and prussian blue. 


colloids,’ such as silver, ferric 
Colloidal ferric hydroxide is precipitated com- 
pletely by pressures as low as 100 atmospheres 
whereas molybdenum blue is only slightly pre- 


cipitated at a pressure of 17,000 atmospheres. 
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The electrode potential! of the hydrogen elec- 
trode was investigated under pressures up to 
thirty thousand atmospheres and the e.m.f. of 
the Weston Standard Cell* up to twelve thousand 
atmospheres. The effect of pressure on both of 
these is of a rather low order of magnetism. The 
effect of pressure upon the penetration of water 
and numerous other liquids into glass and metal 
surfaces” has been investigated with some very 
interesting results. Water was found to penetrate 
in considerable quantities to a depth of several 
millimeters into glass in only a few minutes 
time. The penetration of alcohol and ether is 
somewhat less, and such liquids as paraffin oils, 
glycerin, etc. penetrate scarcely at all. A similar 
effect, but to a considerably lesser degree, was 
observed for the penetration of liquids into 
metals. The penetration of gases into metals® of 
course is more rapid than that of liquid, and 
represents a problem of considerable industrial 
importance, particularly in the case of pene- 
tration of hydrogen into steel!” whereby the 
tensile strength of the steel is reduced to less 
than half its original value. The normal dissoci- 
ation of steam into hydrogen and oxygen, and 
the subsequent removal of the oxygen by its 
reaction with the metal of the high pressure 
steam lines, particularly under conditions where 
a high superheat is used, presents a problem of 
utmost importance in the high pressure high 
superheat steam installations that are bein: 
used at the present time. It has been shown that 
considerable quantities of hydrogen are con- 
tinually escaping through the walls of the hich 
pressure high superheat steam lines, and this 
subject should be thoroughly investigated to 
determine to what extent this is affecting the 
tensile strength of the metal of the steam lines. 
If this effect is cumulative, as is the case under 
the conditions of many of our experiments, it 
might very well become a problem of utmost 
importance in high superheat steam insulation, 
if not indeed a limiting factor in the industry. 
We have found that this effect is not only a 
function of the temperature and pressure, but 
is also affected by the composition and_ heat 
treatment of the steel. It is, therefore, a com- 
bination high pressure and metallurgical inves- 
tigation. 
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Other fields which have been investigated are: 
compressibility of liquids* and solids, the effect 
of pressure upon the viscosity of liquids, and in 
some cases the change of state, or even poly- 
morphic transition,'~® 

A britf summary of the effect of pressure upon 
the properties of matter will reveal some very 
interesting facts. The following characteristics 
for liquids and solids are all affected by pressure : 
density, volume, index of refraction, electrical 
conductivity, thermal conductivity, magnetic 
permeability, dielectric constant, optical rota- 
tion, chemical reactivity, solubility, phospho- 
rescence, fluorescence, physical strength, specific 
heat, latent heat, permeability to gases and 
liquids and viscosity. If we consider all but the 
last two items on this list, we will see that the 
effect is anywhere from a small percent up to 
certainly less than a factor of ten. But now let 
us consider the last two items. The permeability 
to gases and liquids may be of the order of mag- 
nitude of many thousand-fold, and the coefficient 
of viscosity of many lubricating oils is increased 
by as much as two hundred million times their 
value at atmospheric pressure. It is, therefore, 





Fic. 2. Set-up used for absorption spectra measurements 
of materials under pressures of as high as 580,000 pounds 
per square inch. 


our belief that the permeability to gases and 
liquids and viscosity are the two characteristics 
that will represent the most fruitful problems 
for investigation and certainly are the two that 
are most likely to produce noticeable effects in 
industrial processes involving high pressures. 
Investigation of the permeability of gases and 
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liquids has already been discussed and the vis- 
cosity investigation of lubricating oils under 
pressure probably represents the most fruitful 
investigations to be carried out in this field. 





Fic. 3. Assembly used for producing a pressure of 
1,500,000 pounds per square inch. 


Anyone familiar with tests on lubricating oils 
is aware of the fact that a small quantity of 
suspended solids in the oil is very likely to 
produce a bearing failure even though these par- 
ticles may be softer than the two metals being 
lubricated. This effect is in part, if not largely, 
the result of the solid materials which may not 
be so soft at the higher pressures, so completely 
wiping off the lubricant as to permit the two 
metal surfaces to come in contact and start the 
failure. It is the opinion of the author that many 
failures under conditions where an extreme 
pressure lubricant is required are caused by a 
similar phenomena, but that it has not been 
generally recognized because an examination of 
the oil before failure shows no such foreign par- 
ticles, and after the failure only such particles as 
have been removed from the metal surface during 
failure. One very likely source of such suspended 
particles has in general been overlooked, namely 
those which are liquids at normal pressures but 
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which separate as solids under the extreme 
pressure. In some cases, these solids may be of 
metallic hardness. Water in the form of an 
emulsion, or even the normal small amount of 
moisture in the oil is converted to a solid at one- 
third to one-fourth the pressure frequently 
existing in an oil where an extreme pressure 
lubricant is required. The solid thus formed js 
known as ice VI and is very much harder than 
ordinary ice. Under these pressures it has a 
melting point about as high as the normal boiling 
point of water. 

Some of the normal constituents of certain oils 
may go through a similar transition, producing a 
substance hard enough to wipe the lubricant 
from the metal surfaces, or even produce some 
abrasion action itself. Another phenomenon 
which has frequently been observed in our work 
at ultra pressures and which no doubt has its 
effect not only upon the wear, but upon sub- 
sequent failure, is the actual penetration of water 
or other compressed materials into the surface 
of the steel. The escape of such substances when 
the pressure is released causes a breakdown of 
the metal surface.’ If a glass rod in contact with 
water is subjected to such pressures and the 
pressure released immediately, the rod will not 
be affected, but if the pressure is maintained for 
five minutes, so as to allow time for the water to 
penetrate the surface of the glass, and then the 
pressure is released, the rod will be broken up 
into disks; whereas, if the pressure is maintained 
for twenty minutes permitting the water to 
penetrate deeply into the glass, then released 
rapidly, the glass will be shattered into small 
pieces. Similar effects, although to a much lesser 
degree, are known for steel and other metals. 

\With the increasing realization that pressures 
of as high as four hundred thousand pounds per 
square inch exist between two metal surfaces 
which are supposed to be kept separated by 
means of a lubricant, has led to an entirely 
different conception of the problem. Extreme 
pressure lubrication is only in its infancy, and a 
proper method of approach to a fundamental 
understanding of the problem is through research 
at extreme pressures. 

Suppose we put a lubricating oil under such a 
pressure, and study its physical properties. We 
find that it is no longer a fluid in the normal 
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sense of the term, but that it has in some cases conception that a lubricant is a fluid without 
become quite hard. It is indeed fortunate that harmful properties serving as a carrier for a very 
this is true, for if such were not the case, it would small quantity of material which produces a 
be impossible to maintain a film of lubricant lubricating action. The proper hardness curves 
between many of the operating parts of a modern can be theoretically determined from the condi- 
automobile. tion under which the oil is intended to operate, 
The viscosity of a liquid at normal atmospheric — taking into account the hardness and curvature 
pressure is no indication of what its viscosity or of the two surfaces, the rate of motion of the 
hardness may be at greatly elevated pressures. contact point with respect to these two surfaces, 
Two liquids having almost identical physical the rate of motion of one surface with respect to 
properties at One pressure may vary by several the other, and other conditions affecting the 
thousand-fold at another pressure well within — plastic flow of the lubricant. With these speci- 
the range of those normally encountered in fications available, it will be possible for the 
extreme pressure lubrication problems. The lubrication engineer to design a lubricant that 
effect of pressure upon the physical properties of | will have the physical characteristics necessary 
various liquids differs so widely as to make it for it to serve as a lubricant in the true sense of 
possible to produce a mixture of almost any — the word, and which will prevent abrasion of the 
desired viscosity or hardness at almost any two metal surfaces. 
specified pressure. It was the development of a new steel that 
It is the opinion of the author that the day is permitted the production of an automobile 
not far distant when the specifications of a requiring an EP lubricant. One has but to 
lubricant will be largely based upon the pressure consider the present trends in the development 
hardness curves for not only the mixture as a_ of alloy steels to realize that the time is already 
whole, but upon the individual constituents. at hand when an E P lubricant of a still higher 
This is particularly true in view of the modern — order of magnitude is required. 
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Application of the Lagrangian Equations 
to Electrical Circuits 


By D. A. WELLS 


University of Cincinnati, Cincinnati, Ohio 


T has been known 

for many years that 
the Lagrangian equa- 
tions of analytical 
dynamics can be ap- 
plied to electrical cir- 
cuits as a means of 
arriving at the differ- 
ential equations of cur- 
rent flow and mechan- 
ical motion of moving 
parts of the circuit. However, this method has 





never come into general use even though it is 
quite easy to apply and, for certain types of 
problems to be illustrated later, is a remarkably 
powerful tool. There seem to be two reasons 
why the Lagrangian method has not come to 
First, simple 
circuits with no moving parts can be handled 


be used more generally. most 
more quickly by a direct application of Kirch- 
hoff’s laws and it is a fact that most circuits are 
of this type. Secondly, the textbooks of applied 
electricity almost invariably omit a discussion 
of this method. Hence the electrical engineer has, 
in most cases, not had an opportunity to become 
familiar with the use of this powerful tool. 

The purpose of this paper is, therefore, to 
discuss briefly a convenient form of the La- 
grangian equations applicable to electrical cir- 
cuits and to illustrate their use with a number of 
specific examples. A brief review of the applica- 
tion of Lagrange’s equations to mechanical 
systems will first be given. The derivation of 
these equations will not be repeated here since 
it may be found in any 
analytical dynamics. 


advanced text on 
If the generalized coordinates necessary to 
the 
mechanical system are denoted by, 41, qe, -- 


holonomic 


, Qn: 
the Lagrangian equations of motion of the 


represent configuration of a 
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system may be written as 


d / OL aL 
(ter 

dt\ dq, 0g. 
where the symbols have the following meaning: 
gy represents any one of the generalized coordi- 
nates, q, is the first time derivative of g,, and L, 
the “Lagrangian the 
kinetic energy 7 minus the potential energy V 
of the system. That is L=7—JV. F, represents a 
generalized dissipative force such that F,ég, 
gives the energy dissipated by the system when 
g, is changed to q,+46q,, all other coordinates 
remaining unchanged. Since (1) is true for any 
one of the coordinates, g:, ge, * * 


so-called function,” is 


*, Jn, it represents 
a set of m second order differential equations 
which when solved give the complete motion of 
the system. 

As an illustrative example, consider the motion 
of the double pendulum shown in Fig. 1. Let it 
be assumed that the motion is confined to a 
plane, that the lengths 7; and 72 are fixed and 
that a damping vane is attached to m,. It is 
that 
coordinates are required to determine completely 


seen under these conditions only two 
the positions of the masses wz, and mp. Of the 
many choices of coordinates that could be made, 
the angles @ and ¢ are perhaps the most con- 
venient. We shall now set up the Lagrangian 
function ZL and by an application of (1) write 
out in explicit form the two differential equations 
of motion of the pendulum. 
The kinetic energy of mm, is seen to be 


T = (m, 2)(r6)?. 


The kinetic energy of mz can most easily be 
found by first writing down its position in the 
From the diagram it 


x and y coordinates. 
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Fic. 1. A double pendulum for illustrating the application 
of Lagrange’s equations to particle dynamics. 


follows that 
x=r, sin 0+r2 sin ¢+%Xo, 


y=r,(1—cos 0)+r2(1—cos 9) +0, 


where x» and yo are the coordinates of mz in its 
rest position. Hence, the x and y components of 
velocity are 

i=7,6 cos 0+17re¢ cos ¢, 

y=7,6 sin 0+re¢ sin ¢. 


The kinetic energy of mz is, therefore, 


T= (me/2)(4?+9") 
= (me/2)[71°6? + re? ge? + 2r robe cos (¢ —8) | 


ad 
and the total kinetic energy is given by 
T=7T,+72=3[(mi4+mye)r)°e? 


+ more? ge? + 2morire¢ cos (¢—9) | 


J 


The potential energy of the system, referred to 
the lowest possible position of mo, is easily seen 
to be 


V=(m,+me)gr,(1—cos 0) +mogro(1—cos ¢). 


Assuming that the viscous drag on the vane 
attached to m, is proportional to its linear 
velocity (the proportionality factor being k), 
the change in energy of the system by a small 
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displacement, 7:66, of m, is given by 
6& = — (kr16)(7,60) = F140. 


Hence the generalized dissipative force F, corre- 
sponding to the coordinate @ is 


F, = kr,°6. 


Since the damping of ms. is assumed to be 
negligible, F,.=0. The Lagrangian function is 
now seen to be 


L=3[(mi4+mz)r°6?+mor2¢? 
+2moer rb cos (¢—4) } 
— (m,+mye)gr;(1—cos 6) — megre(1—cos ¢). 
From which, 
(OL / 00) = (m,+mz2)r°6+ moerire¢ cos (y—8), 
(d/dt)(0L/ 00) = (m,+mz2)r,26+morireg cos (¢—8) 
+ mor \ro¢e(6— ¢) sin (¢—8), 
(0L/00) = mor reb¢ sin (g—4) 
—(m,+me)gr; sin 0. 
Therefore, by (1), the first equation of motion is 
(m,+m2)r 126+ mor reg cos (g—4A) 
+ mor ree(6— ) sin (¢—0) —moerrob¢ sin (g—8) 
+ (m,+mz2)gr; sin = —kr,6. 


Likewise the second equation of motion is found 
to be 


More? e+ mor rb cos (y—8#) 
+ mor \r26(6— ¢) sin (¢—8) + mor rb¢ sin (¢—A) 
+moegre sin g=0. 


When ZL and the dissipative forces F, have 
been properly determined, Lagrange’s equations 
give the differential equations of motion of any 
mechanical system in any system of coordinates. 
(It should be noted that when L and F, have 
been determined the labor of determining the 
equations of motion is over. The final step or 
the application of (1) amounts only to a bit of 
simple differentiation.) 

Likewise it is possible to find a Lagrangian 
function and dissipative that 


‘forces’? such 
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Lagrange’s equations will give the differential 
equations of an electrical circuit, however com- 
plex, even if various parts of the circuit are allowed 
to move in any way with respect to each other. 


The Lagrangian Function for an Electrical Circuit 


If (1) is to be used as a means of arriving at 
the differential equations of an electrical circuit, 
it follows that we must find a suitable Lagrangian 
function L and the proper dissipative forces F,. 
This implies a choice of ‘‘coordinates’’ and a 
determination of the ‘‘kinetic’’ and ‘‘potential”’ 
energy of the system. To begin with, it will be 
assumed that no part of the circuit is allowed 
to move. Thus ‘‘mechanical”’ kinetic and poten- 
tial energies will not enter the first part of the 
discussion. 

An explanation of the choice of coordinates 
here made cannot be given in a few words, but 
it will become clear from what follows that the 
charges which, from some initial time f), have 
passed through the various branches of the 
circuit constitute an appropriate set of coordi- 
nates. Hence the charges, represented by Q), Qe, 

-, Q,, correspond in an electrical circuit to 
the coordinates qi, qe, *, Qn, in a mechanical 
system. A fact of great importance in the use of 
charges as generalized coordinates is that there 
are almost invariably relations between certain 
of the charges. That is, not all of the coordinates 
are independent. At any junction in a circuit 
the algebraic sum of the charges flowing to the 
junction must be zero (Kirchhoff'’s law). Hence 
an application of this law gives a set of equations, 
by means of which, certain coordinates (charges) 
can be eliminated. For example, an unbalanced 
Wheatstone bridge circuit contains six branches 
and therefore there are six currents and _ six 
corresponding charges. However, there are three 
independent relations between these charges and 
hence, by elimination, the six can be reduced to 
three. It is a matter of indifference as to which 
three are eliminated. In the language of ana- 
lytical dynamics, the Wheatstone bridge is a 
problem of three ‘degrees of freedom’’—not six 
as might first be supposed. By this method an 
electrical circuit must be reduced to its proper 
number of degrees of freedom. If this is not done, 
Lagrange’s equations, in the form in which we 
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are to use them, will not apply. At the outset of 
work on a circuit each current and each electro- 
motive force must be assigned a direction just as 
if the circuit were to be handled by the usual 
methods. This makes it possible to determine the 
proper number of degrees of freedom and, as 
will be seen later, the potential energy of the 
sources. 

A consideration of the magnetic energy of an 
inductance will show, at once, what corresponds 
to kinetic energy in the electrical circuit. The 
energy of the magnetic field about a coil of self- 
inductance MV, is WJ? 2 where J=(dQ dt) =Q is 
the current in the coil. Since Q is considered as a 
coordinate, Q is a generalized velocity and it is 
seen that ./Q* 2 has the same form as the 
familiar term mv* 2 representing the kinetic 
energy of a particle. 

Consider two coils having self-inductances of 
MJ, and Ms and each a mutual inductance with 
respect to the other of /;.. Denoting the current 
through the first by Q; and through the second 
by Qs, the energy of inductance is given by 


MQ? M.Q.? 
= + 


2 


= - 


g + MW y2QiQ>. 


But this is again a quadratic function of the 
generalized velocities and corresponds closely in 
form to the expression for the kinetic energy of 
the double pendulum. A further consideration of 
the energy of inductance shows that it can 
always be expressed as a quadratic function of 
the currents or generalized velocities. Hence, 
energy of inductance is to be regarded as kinetic 
energy. 

The potential energy of a circuit is made up 
of two parts, the energy of the sources (batteries 
etc.), and the energy stored in condensers that 
may be connected in the circuit. A source will 
always deliver, or tend to deliver, charge to 
the circuit. If a charge Q is delivered by a source 
having a constant electromotive force EF, the 
source has less energy to the extent of EQ. 
Regarding the source as having potential energy 
and referring it to the “‘point’’ Q=0, it is seen 
that this potential energy must be written as 
— EQ. This amounts to saying that in order to 
reduce a certain delivered charge to zero, work 
must be done on the source. Of course, if a charge 
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is regarded as flowing against the electromotive 
force of the source, the potential energy EQ must 
be considered positive. That the potential energy 
of a source with constant electromotive force is 
exactly analogous to the potential energy of a 
particle of mass m acted upon by gravity, is 
easily seen. If y, Fig. 2, is considered positive 
down, the potential energy of m is —mgy, mg 
corresponding to & and y to Q. It is important 
to note that even though the electromotive force 
of a source may vary with time, as for example 
E=E,) sin wt, —EQ is still to be regarded as a 
potential energy. This is true for the reason that, 
in the derivation of Lagrange’s equations, time 
is not allowed to vary as any one of the coordi- 
nates g, is changed to q,+6q,.! 

The energy of a charged condenser is Q*? 2C, 
where C is the capacity. Since this expression 
has exactly the same form as the expression for 
the potential energy of a stretched spring, it is 
immediately suggested that the energy of a 
charged condenser be regarded as _ potential 
energy. It is clear that this energy must be 
taken as positive since, when the charge is 
reduced to zero, energy is expended by the 
condenser—not on it. Hence the potential energy 
of a circuit is the algebraic sum of the energies 
of the sources plus the sum of the energies of the 
condensers. 

It only remains to determine what is to be 
regarded as the generalized dissipative forces. 
Neglecting radiation, hysteresis, and eddy cur- 
rents, the electrical energy of a circuit can only 
be dissipated as heat in the resistances. Of 
course, electrical energy can be changed to 
mechanical energy by allowing parts of the 
circuit to move but, as previously mentioned, 
we are for the moment assuming that each 
part of the electrical system is fixed. The rate 
of change of energy of a circuit containing only 
one resistance R, due to the passage of current 
through this resistance, is given by 


(de dt) = —RQ?. 


But FQ, where F is the generalized force and Q 
the current through the resistance or generalized 
velocity, must also represent the rate of change 


‘See Whittaker, Analytical Dynamics, p. 36. 
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of energy of the system. Hence, 
FQ= — RQ. 


Thus the generalized dissipative force is deter- 
mined by 
F= —RQ. 


A complex circuit requires as many differential 
equations for its treatment as there are degrees 


asi 





+O 


mg 





Fic, 2. Simple diagram to illustrate the analogy between 
gravitational potential energy and the potential energy 
of an electrical source. 


of freedom (independent currents or charges). 
Each equation involves a generalized dissipative 
force and as a general rule each force involves a 
number of currents and resistances. The method 
of arriving at dissipative forces for a complex 
circuit is quite easy. It is best shown by the 
examples that follow. 

We, therefore, write the Lagrangian equations 
applicable to a stationary electrical circuit as 


d sol aL 
Fn ee | 
dt aq). dQ. 


where Q, represents any one of the nm charges 


Nm 


necessary for the determination of all other 
charges of the circuit, Q. the corresponding 
current, F, the corresponding electrical dissipa- 
tive force and where, in the Lagrangian function 
L=T—TI, the kinetic energy T and potential 
energy | are determined as outlined above. 

In order to make clear just how the L function 
is set up and how the circuit equations are ob- 
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Fic. 3. Simple series circuit. 


tained from the Lagrangian equations, two 
specific examples will now be given. Consider 
first the simple series circuit of Fig. 3. The 
kinetic energy T=MQ?/2 and the potential 
energy V = —EQ+(Q*/2C. Hence the Lagrangian 
function appropriate to this circuit is 


MQ? Q? 
L=T-—V=——+EQ-—. 
2 2C 
The dissipative force F= —RQ. Therefore, since 
d saL - aL O 
(=) =MQ and —=E-—-, 
dt\aQ dQ C 


an application of (2) gives 
MO-—E+0/C=—RQ, 


which is the familiar differential equation of this 
circuit. 

Consider next the simple circuit of Fig. 4. 
We shall assume that the inductance coils each 
have a mutual inductance of M2 for the other. 
The charges to be considered are Q, Q;, and Qs. 
But since Q= Q:1+(Qz, one of these charges, say Q, 
can be eliminated. Hence the kinetic and poten- 
tial energies are given by 


— MQ M2Q.? _ 
T =——- ++ M2002, 
2 2 
—_ rs 
V=—E,(Q:+Q2) +£202+—+—. 
2C, 2C2 


Since this is a problem of two degrees of freedom, 
the coordinates being Q; and Qe, there are two 
corresponding generalized dissipative forces F 
and F, which can easily be found as follows. 
Suppose the charge Q is thought of as increasing 
slightly. The change in energy of the system is 
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given by 
66 = — R,Q16Q1 — R2Q26Q2— Rs(Q: +Q2) (50: +40,) 
= —[RiQ:+R3(Q:+Q2) ]60; 
—[RQ2+Rs(Qi+Q2) 40x, 
Hence, for a slight change of Q, only, 


66:= —[RiQi+Rs(Qi+Qs) JQ). 


Setting this equal to the generalized force F, 
times the displacement 6Q, we get 


Fi = —[RiQi+R3(Q:+Q.)]. 
Likewise F.= —[RoQ2+R3(Qi:+Q.) ]. 


d so0L OL 
Thus, from — ( ~) ———_=F 
d dQ dQ; 
d s aL OL 
and -(—) ——— = F,, 
dt OQ» dQ. 


we obtain the two differential equations of 
“‘motion”’ of the circuit 


M,0,+ Mi202—E,+(01/C)) 

= — R,Qi—R3(Qi+Q:), 
M2Q2+ M 120: — E:+E2+(Q2/Ce2) 

= —R.Q2—R;3(Q:+Q). 


These differential equations and the one for 
the first example are, of course, just the equations 
usually obtained by an application of Kirchhoft’s 
laws. It is seen that the Lagrangian method has 
little if any advantage when applied to the simple 
type of circuit illustrated thus far. The real 
power and usefulness of the method becomes 
evident when applied to circuits where certain 
parts are allowed to move. Following a brief 
discussion of the form taken by the Lagrangian 
function for such a circuit, two further examples 
will be given. 

If parts of an electrical circuit are in motion, 
the total kinetic energy of the system must 
include the kinetic energy of motion of the parts 
as well as the energy of inductance. Likewise the 
expression for the potential energy must include 
terms for the potential energy that the circuit, 
treated as a mechanical system, may have as 
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well as the potential energy of the sources and 
energy of the condensers. In other words, the 
Lagrangian function must contain the kinetic and 
potential energies of a mechanical system as well as 
the electrical circuit. This means that the La- 
grangian contains two groups of generalized 
coordinates, the charges Q:, Qe, ---, Q,, and the 
usual displacement coordinates qi, g2, --*, Qm- 
The number of differential equations now ob- 
tained and the number of generalized dissipative 
forces required will be n+m. If there are no 
forces of friction, viscous drag, etc., acting on 
the moving parts the generalized forces corre- 
sponding to qi, d2, ***, Ym, are each zero. 

As a simple example of an electrical system 
containing moving parts consider the circuit 
shown in Fig. 5. The upper plate of the condenser 
C, having a mass of m, is attached to the coil 
spring and is free to move in a vertical direction 
under the action of the spring and electrical 
forces. Gravity is not considered. All other parts 
of the circuit are fixed. The unusual features of 
this circuit are the variable capacity C, the 
kinetic energy of the moving plate, and the po- 
tential energy of the spring. 

Call s the distance between condenser plates 
when the spring is unstretched and x the dis- 
placement of the plate from the position it 
occupies when the spring is unstretched. The 
potential energy of the spring is thus kx*/2, 
where k is the usual spring constant. The capacity 
of the condenser is given by 


A a 


C= —= ; 
4n(st+x) s+x 





where A is the area of one plate and the dielectric 
constant is taken as unity. The potential energy 
of the condenser is, therefore, (Q?/2a)(s+x) and 
the kinetic energy of the moving plate is mz?/2. 
Thus the Lagrangian for this system is seen to be 


9 


MQ?) mz? Q? : 
L=——+—+EQ-——(s+x) -—-—. 
2 2 2a 2 


This is a problem of two degrees of freedom, the 
two generalized coordinates being Q and x. 
From the general Lagrangian equation the two 


differential equations of the circuit follow at 
once : 
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MOQ+(Q/a)(s+x)-E=RQ, 
mé — Q?/2a+kx=0. 


A solution of these two simultaneous differential 
equations gives the mechanical behavior of the 
moving plate as well as the electrical behavior 
of the circuit. If a mechanical dissipative force 
is allowed to act on the moving plate it can 
easily be taken into account. For example, if a 
viscous drag (given by f=—bz) is acting on 
the moving plate, the second equation becomes 


mé—Q?/2a+kx = — be. 


In the more general problem any part of the 
circuit may be allowed to move. As a final ex- 
ample and one which enables us to write down 
the Lagrangian function in a very general way, 
consider the circuit of Fig. 6. Due to the length 
of details that would be involved several steps in 
the work are left in implicit form. It is believed, 
however, that no details are omitted that would 
cause great difficulty in handling an actual 
problem. 

There are four charges to be considered but 
they are not all independent. It is seen from the 




















Fic. 4. Circuit to illustrate the method of obtaining 
‘dissipative forces.” 


diagram, assuming that the charges are flowing 
in the directions indicated, that Q;-—Qs=Q2—(Q). 
Hence any one of the charges can be eliminated 
and the problem, so far as the electrical part is 
concerned, is one of three degrees of freedom. 
In the relations that follow we will eliminate Q,. 

Assuming that the coils each have the self- 
inductance indicated on the diagram and that 
each has a mutual inductance with respect to 


317 





r 
™ 
TOON 


M 


Fic. 5. Circuit containing a moving plate condenser. 


every other, the energy of inductance or the 
energy corresponding to kinetic energy is given 
by 


Teroet. = 3(M11Q12? + M2Q1Q2+ M1301Qs) 
+3(Mo1QoQ1 + Mo2Q2? + Mo3Q2Qs) 
+3 (M51Q3Q1 + M52Q:Q2+ Ms3Q;°), 
where, for example, 1/2; is the mutual inductance 


of the second and third coils. This can be written 
in compact form as 


3 
T ctect. = 43>. M,.Q,Q,. 


The potential energy of the condensers is 


QO? Q:? Q% Qe 
+—+—+—. 


Q271 1\ Qky1 1 
Mele) 
2AC, Cs 2A\C2 Cs 
Q3;771 = 1 2:03 — Q203— 2102 
(142). oa 
C; 
The potential energy of the sources is 
Sr = — E,Q,— E2Q2— E3Q3+ EsQ. 
Eliminating Q; we have 
Ee=Qi(Ey— E;) —Qo(E2+ Es) +Q3(Es— Es). 
The electrical dissipative forces are found as 
follows. Suppose that each of the charges is 


given an infinitesimal change, subject to the 


conditions that Q;=Q;+Q:1—Qe, 6Q0;=603;+ 40, 
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—6Q2, and Q;=Q;+Q:—Q2. The change jp 
energy of the system due to heat produced jp 
the resistances by these small changes in charge 
is given by 


66 = —[RiQ1601 + R2Q2602+R3Q36Q0s +R .Q,6Q,] 
Eliminating 6Q, and Q; we get 
66 = — {60,(RiQi + Rs(Qi+Qs—Qz) ] 
+ 602 R2Q2— Ri(Qi+Q;—Q:) ] 
+60s[RsQ:+Ri(QitQs—Q2)}}, 


Therefore, for a displacement of 6Q; only, we 
have 


66, = F,6Q,= — 62:[RiQi+Rs(Qit+Qs—Q.)]. 


Hence the generalized dissipative force corre- 
sponding to Q; is 


Fi = —R,Q:—Rs(Qi:+Q;—Q.). 
Likewise, 

F2= — RoQo+Rs(Q3+Qi—Qo) 

F;3= —R3Q3—Rs(Q;+Q:—Q,). 


If the inductances are free to move and if, 
perhaps, certain condenser plates are allowed to 
move, then the circuit regarded as a mechanical 
system possesses kinetic and potential energy. 
As previously mentioned, this energy must be 
included in the Lagrangian function. In order to 
represent this energy we must, of course, intro- 
duce a set of coordinates q;, ge, «++, Gm, sufficient 
in number to represent completely the con- 
figuration of the system. (If one coil only were 
allowed to move without rotation in the plane 
of the paper, then only two coordinates, say x 
and y, would be required. If, in addition to this, 
one condenser plate were allowed to move, 
three displacement coordinates would be re- 
quired.) 

A general expression for the kinetic energy 
of a mechanical system referred to fixed axes is 

m 


a — — 
1 mech. — 2 da rs7rQs; 
rs 


where ¢ represents a generalized velocity corre- 
sponding to the generalized coordinate qg and 
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where the coefficients a,,, containing the masses 
and in some cases the moment of inertia of the 
parts, are functions of the displacement coordi- 
nates. The indices of the double sum each run 
from 1 to m, m being the number of displacement 
coordinates necessary to represent the configura- 
tion of the system. Referring to the double 
pendulum, m=2, qi=9, ga=¢, di=(m,+m2)r;’, 
dog= Mors”, ANd Aj2=A21= MePiPr2 COS (Y—84). 
The potential energy of the system, treated as 
a mechanical system, will depend entirely on the 
particular circumstances and can only be written 
explicitly when these circumstances are given. 
For example, if a set of springs is holding the 
coils and condenser plates and if, perhaps, 
certain parts are allowed to move under gravity, 
the potential energy can be written in explicit 
form. If the circuit is on a horizontal table with 
no springs or other conservative forces acting on 
the parts, the potential energy is, of course, zero. 
We shall denote the potential energy of the 
circuit, treated as a mechanical system, as V’. 
In like manner the dissipative forces corre- 
sponding to the various coordinates qi, qo, 
‘, Gm, can only be stated explicitly when 
definite conditions are known. The dissipative 
é, M 
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Fic. 6. Circuit to illustrate the more general form of 
Lagrange’s equations. All condenser plates and all in- 
ductances may be allowed to move. 


force acting on m, of the double pendulum is a 
typical example. Generalized dissipative forces 
corresponding to the displacement coordinates 
will be denoted by f1, fs, «++, fm. In any actual 


problem any one or all of these may be zero. 
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We are now in a position to write down the 
Lagrangian function for the circuit. Including 
both electrical and mechanical energies, it is 
clear that 


3 m 
L=31>M,.0,0.+3D.0,.6-4.—-Ec—Er—V’. 


The Lagrangian equations for the system are, 
therefore, 


d fo0L OL 
- ( : ) —_ = F, ’ (3) 
dt\ dQ, dQ, 
d f/dL OL 

( , ) -—=f,, (4) 
dt\ 0q, 0g; 


where, for the particular circuit in hand, r in 
(3) extends from 1 to 3 and in (4), from 1 to m. 
When the operations of (3) and (4) are carried 
out there result 3+-m second order differential 
equations. A solution of these simultaneous 
equations gives the complete electrical behavior 
as well as mechanical motion of the circuit. 
It is important to note that when inductance 
coils and condenser plates are allowed to move 
the expressions for mutual inductance and 
capacity become functions of the coordinates 
Gi, G2, ***, Gm. This fact must be taken into 
consideration when performing the operations 
of (3) and (4). 

In the more general case where a circuit is 
composed of any number of inductances, re- 
sistances, and condensers, any one or all of 
which are allowed to move, it is seen that the 
appropriate Lagrangian function is 


tw 


) 


n m 
L=31DM,,0,0.4+3D0n464.+ DL EO-1E—-V’, 
where it is assumed that, before L is made use 
of (see the previous example), 
charges have reduced to a_ necessary 
minimum by Kirchhoff's law and also the 
displacement coordinates have likewise been 
reduced to the necessary minimum by the con- 


the number of 
been 


straints on the system. 

It is clear from the first two examples given, 
that the Lagrangian method of arriving at the 
differential equations of a simple stationary 
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electrical circuit offers no particular advantage 
over the usual Kirchhoff’s law method. How- 
ever, for a combination electrical-mechanical 
system, particularly where it may be advan- 
tageous to use coordinates for the moving parts 
other than the simple rectangular type (La- 
grange’s equations apply to any system of 
coordinates), the Lagrangian method reduces 
what otherwise might be very difficult to a 
more or less routine task. When the Lagrangian 


function L has been properly determined, terms 
in the final differential equations are auto. 
matically given the proper algebraic signs. Like- 
wise all interrelations of electrical and mechanica] 
quantities are automatically taken care of. 

Since electrical engineering and its related 
branches are replete with combination electrica]- 
mechanical problems, it is believed that this 
method has a wide field of usefulness in modern 
engineering. 





New Appointments and Promotions 


Dr. H. R. Crane has been promoted to an assistant 
professorship in the department of physics at the Uni- 
versity of Michigan. 


* 


Professor John A. Wheeler of the University of North 
Carolina has been appointed assistant professor of physics 
at Princeton University. 


* 


Dr. Millard F. Manning of the University of Toledo 
has accepted an appointment as assistant professor of 
physics at the University of Pittsburgh. Professor Manning 
will teach primarily courses in the structure of metals. 


* 


At the Massachusetts Institute of Technology professors 
Robley D. Evans and Joseph C. Boyce have been pro- 
moted from assistant professors to associate professors; 
Drs. Walter Albertson, L. C. Van Atta, and C. M. Van 
Atta have been promoted from research associates to 
assistant professors; Dr. G. G. Harvey of the College of 
the City of New York has been appointed assistant pro- 
fessor; and Mr. J. H. Buck of the University of Rochester 
has been appointed instructor. 


7 


L. P. Delsasso has been granted a leave of absence for 
one year to serve as exhibits expert in the physical sciences 
for the University of California exhibit in the Golden 
Gate International Exposition. 
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Professor Maurice Ewing has been awarded a Guggen 
heim fellowship and will be on leave of absence from 
Lehigh University during 1938-1939. 


* 


The Burgess Battery Company has elected Mr. Dan W. 
Hirtle president and Dr. C. F. Burgess chairman of the 
board of directors at its annual meeting held March 26, 
1938. Prior to his association with Burgess, Mr. Hirtle 
was an assistant superintendent of the American Rolling 
Mills Company, Middletown, Ohio. 


* 


Dr. Lee A. DuBridge, Harris professor of physics and 
chairman of the department, has been appointed dean of 
the faculty of arts and sciences of the University of 
Rochester. He succeeds Dr. Leonard Carmichael, who has 
been elected president of Tufts College. 


* 


Summer Plans 


Professor I. I. Rabi of Columbia University will be a 
member of the summer quarter staff at Stanford Uni- 
versity. 


* 
Sigma Xi Lecture 


On Saturday, April 2, Dr. W. F. G. Swann, director of 
the Bartol Research Foundation, gave the annual Sigma 
Xi lecture at Carnegie Institute of Technology. 
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Research Fellowship Appointments Announced 


The first five Westinghouse Research Fellows have been 
selected for work in fundamental science at the Westing- 
house Research Laboratories in East Pittsburgh. The 
successful candidates chosen from a group of fifty are: 


Robert O. Haxby, Univer- 
sity of Minnesota. Mr. 
Haxby, whose home is in 
Minneapolis, is twenty-six, 
has received the bachelor’s 
and master’s degrees in 
electrical engineering from 
the University of Minne- 
sota. He has been a grad- 
uate student in _ physics 
there for the last four 
years and has participated 
actively in the nuclear 
physics research there. He 
expects to receive his doctor’s degree in physics next June. 
He will join the group who will carry on research in nuclear 
physics, using the large high voltage generator now being 
built in East Pittsburgh. 





John A. Hipple, Prince- 
ton University. Mr. Hipple 
is twenty-seven. His home 
is in Lancaster, Pennsyl- 
vania, where he received a 
bachelor’s degree in physics 
from FrankJin and Marshall 
College. He has been a 
graduate student in physics 
at Princeton during the 
past four years, where he 
has worked with Professor 
Walker Bleakney and where 





he expects to receive the de- 

gree of doctor of philosophy in June. He will design and 
construct a large mass spectrograph of high resolving power 
for studies in molecular structure in the Westinghouse 
Laboratories. 


W. E. Stephens, Cali- 
fornia Institute of Tech- 
nology. Stephens is twenty- 
six, and expects to receive 
the doctor’s degree next 
June for researches in nu- 
clear physics conducted in 
the Kellogg Radiation Lab- 
oratory in Pasadena. While 
in Pasadena, he has collab- 
orated with Professor C. 





C. Lauritsen. He is a native 
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of St. Louis, Missouri, where he received a_ bachelor’s 
degree in 1932 and a master’s degree in 1934 from Wash- 
ington University. At Washington University he did 
researches on focusing of ion and electron beams in col- 
laboration with Professor A. L. Hughes. In Pasadena, 
Stephens has studied the production of neutrons by bom- 
barding boron and nitrogen with high energy deuterons. 
He has also worked on the development of a new type of 
high voltage discharge tube. In East Pittsburgh, he will 
continue his researches on nuclear physics with the large 
generator at the Westinghouse Laboratories. 


Sidney Siegel, Columbia 
University. Dr. Siegel is 
twenty-six and received his 
doctor’s degree two years 
ago for researches on elastic 
and magnetic properties of 
metals carried out under the 
direction of Professor S. L. 
Quimby. Dr. Siegel is a na- 
tive of New York Citv, and 
received his bachelor’s de- 
gree in 1932 from Columbia 
College. Dr. Siegel will con- 
tinue his studies on funda- 





mental properties of metals. 
He plans to develop methods of growing single crystals of 
pure metals, and of alloys which show order-disorder tran- 
sitions in order to study how the elastic properties of alloys 
depend on the arrangement of the atoms in the crystal 
lattice. 


W. E. Shoupp, Univer- 
sity of Illinois. Dr. Shoupp 
is at present an instructor 
in physics at the University 
of Illinois, where he is 
actively engaged in studies 
on nuclear physics, with 
the Illinois cyclotron. Dr. 
Shoupp is thirty, and his 
home is in Troy, Ohio. He 
received his bachelor’s de- 
gree from Miami Univer- 
sity, Oxford, Ohio, and 
went to Illinois for graduate 
study in 1931. In addition 





to his work in nuclear physics, he has made several 
contributions to extreme ultraviolet spectroscopy as a 
collaborator of Professor P. G. Kruger. At the Westing- 
house laboratories, he will continue studies in nuclear 
physics, working with Dr. W. H. Wells on the high voltage 
generator. 


The five appointments are the first to be made under a 
plan announced last December, whereby the Westinghouse 
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Electric and Manufacturing Company plans to support 
the work of ten young physicists on fundamental studies 
broadly related to the electrical industry. The fellows are 
appointed for one year and are eligible for one reappoint- 
ment. They will work under the general supervision of Dr. 
E. U. Condon, Associate Director of the Laboratories. 

“We have been extremely gratified by the response to 
the fellowship plan which is represented by the splendid 
group of applications received,” said Dr. L. W. Chubb, 
Director of the Research Laboratories, in discussing the 
appointments of the first group of fellows. ‘‘This made the 
task of selecting the men a difficult one. We regret only 
that it is not possible to offer appointments to a larger 
number of the applicants.” 


* 
The Electrochemical Society 


The Electrochemical Society held its Spring meeting 
at Savannah, Georgia, April 27 to April 30, 1938. The 
Society met jointly with paper experts to review and 
discuss old and new problems with which that industry is 
confronted. The symposia were scheduled on industrial 
applications of electro-osmotic phenomena and on the 
electro-deposition of metals, particularly manganese, 
chromium, zinc, copper, nickel and alloys of copper and 
of nickel. A fourth session was devoted to principles of 
electrochemistry including papers on the dropping mercury 
electrode, pH determination, new differential dyna- 
mometer, motor-electrolytic potentials, film potentials, 
and over-voltage. 


* 
Institute of Radio Engineers 


The Institute of Radio Engineers cordially invites all 
interested physicists and engineers to attend its Annual 
Convention at the Hotel Pennsylvania, New York, June 
16th-18th. At the technical sessions, fifty papers are 
scheduled for presentation and discussion on vacuum tubes, 
communication circuits, television, and other subjects of 
value to workers outside the radio field. 

In addition, a comprehensive display of laboratory in- 
struments, communication-engineering equipment, vacuum 
tubes, and materials will be open for inspection throughout 
the convention. Thirty leading manufacturers are partici- 
pating in this exhibition. 

The June issue of the Proceedings of the I. R. E. will 
carry the complete program with summaries of the tech- 
nical papers or a list of papers to be presented may be 
obtained from the I. R. E. Convention Secretary. There 
will be no registration fee or other charge. 


* 


The Institute of Metals, London, has made the first 
award of its new medal to Sir. William Bragg, president 
of the Royal Society. The presentation was made on 
March 8. The medal—in platinum 
standing services to nonferrous metallurgy. 


is given for out- 
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New Research Division 


The formation of a division of engineering and research 
by Crane Co., Chicago, with L. W. Wallace, head of 
engineering research for the Association of American 
Railroads, as director has been announced by Charles BR 
Nolte, president of Crane Co. The new Crane division 
will comprise the existing division of research and de. 
velopment and the product engineering department of 
the company. 

“The new division has been formed to coordinate ajj 
engineering activities of the company,” said Mr. Nolte, 
“and to further its progress in diversified fields. The 
importance and growth of this technical work have proved 
the advisability of maintaining a complete and well-staffed 
engineering division, free from all responsibilities except 
those relating to engineering, research, design and experi- 
mental work.” 


* 


Engineering Lectures 


Dr. Reinhold Rudenberg, formerly chief electrical 
engineer of the Siemens-Schuckert Company and honorary 
professor at the Technische Hochschule in Charlottenburg, 
now with the General Electric Company of England, gave 
a series of three lectures at Harvard University on the 
following dates: 

May 11, Technique and economy of electric high power transmission 
over long distances. 

May 12, Oscillations common to synchronous generators and asynchro- 
nous motors in disturbed networks. 


May 13, Straining of machine windings by short circuit forces and 
traveling waves. 


* 


Industrial Research Institute 


A new Industrial Research Institute is being organized 
by the Division of Engineering and Industrial Research 
of the National Research Council. The purpose of this 
institute is to consider general problems which will be of 
interest to a large number of research laboratories such as 
the selection, training, handling, and incentive and bonus 
plans for research personnel. They plan also to consider 
the problem of securing adequately trained industrial 
research workers, the organization, administration and 
operation of budgets of the research laboratories in relation 
to that of the other departments of a company, and the 
coordination of research activities with market studies 
and sales promotion. Recently an executive committee 
was elected with Robert B. Colgate as chairman. Further 
information concerning this institute may be obtained 
from Maurice Holland, who is serving as executive officer. 


* 


Necrology 


We regret to announce that Professor P. A. Ross of 
Stanford University died March 13, 1938. 
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RCA Television Parts 


RCA has recently made available certain specialized 
television parts for use by experimenters within radius of 
television transmitting stations. With these new parts and 
other standard parts already available, it is now possible 
for the amateur experimenter who is equipped with suffi- 
cient technical knowledge, to assemble his own Kinescope 
deflecting circuits for use in experimental television re- 
ceivers. The new parts listed for sale include a deflecting 
voke, two power transformers, a vertical output reactor, a 
vertical oscillation transformer, a horizontal oscillation 
transformer, a horizontal output transformer, two power 
supply capacitors, and a power supply reactor. A folder 
describing these parts and giving wiring diagrams is 
available from RCA Parts distributors or from the RCA 
Parts Division, Camden, New Jersey. 


* 
Photo-Tube Relay 


Built for outdoor service 
in all kinds of weather, this 
is a new photo-tube relay 
unit developed by G-M Lab- 
oratories, Chicago. Complete 
protection from rain and snow 
is given by the heavy water- 
tight case and lens visor. A 
specially designed resistor is 
mounted in the top of the 
case, with fins for dissipating 
heat. The unit is designed 
for 115 volt d.c. operation. 


* 


d’Arsonval Motion Picture 


The French producer, Jean Benoit-Levy, has released a 
moving picture on d’Arsonval’s discoveries. Although 
d’Arsonval is over eighty-four vears of age, he plays his 
own title role. Professor Barail presented the picture to 
the public on Wednesday, April 20, 1938, at the McMillin 
Academic Theater of Columbia University. The per- 
formance was given under the auspices of the Institute of 
Arts and Sciences and the Institut des Etudes Frangaises 
of Columbia University. 


* 


Solving Machine-Tool’s Special Problems 


The Westinghouse Technical Press Service announces a 
squirrel-cage induction motor with an adjustable no-lead 
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speed. Ordinary induction motors come quickly to almost 
synchronous speed, and voltage control or high rotor re- 
sistance do not appreciably affect the light-load speed. A 
novel motor and control has been devised which accom- 
plishes the very unusual result of permitting an induction 
motor to be run at light-load at any desired low speed. The 
scheme consisted essentially of internally loading the motor. 
The motor has a three-phase delta winding, with one phase 
disconnected from the other two, making it essentially an 
open-delta motor. When some lower light-load speed is 
desired the third phase is connected to the line in reverse. 
The current through it is in a generating instead of motor- 
ing direction and tends to buck the motor action of the 
two open-delta phases. The motor slows down as though a 
load had been applied to motor with a high resistance 
rotor. By controlling with a rheostat the current through 
the reversed phase the bucking action and, thereupon, the 
motor speed can be closely controlled. This motor has been 
devised especially for chucking work into a large lathe 
where it is necessary to close in the chuck jaws firmly but 
without shock. 


* 


New Recticharger 
The Raytheon Manufacturing Company announces a 
new automatic charging device which does everything 
necessary to keep a storage battery fully charged except 





BATTERY ALWAYS Fut 


adding the water to make up for evaporation. The control 
unit is designed to adjust the charging current rate in 
proportion to the current drawn from the battery. Over- 
charging of the battery is prevented by this unit and the 
battery will be kept fully charged. The rectifier unit in 
this charger is a copper oxide one. The complete charger 
is mounted in a metal cabinet. 
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Artificial Lightning at World’s Fair 


General Electric’s high voltage laboratory at Pittsfield, 
Massachusetts, where man-made lightning has been 
demonstrated will be transported to the New York 
World’s Fair in 1939. The general public will have an 
opportunity to see in action this artificial phenomenon, 
created for the purpose of duplicating the forces of natural 
lightning. The laboratory will be housed in Steinmetz 
Hall, named after the late Dr. Charles P. Steinmetz, 
credited with being the father of man-made lightning. 
Present plans call for the display of a 10,000,000-volt bolt 
of lightning leaping across a gap of 30 feet with a crash 
greatly simulating the thunder which accompanies natural 





lightning. There will also be the spectacular million- -volt, 
three phase arc, whose ribbons of electrical fire will twist 
their way for 30 feet in the air above the apparatus 
There will also be the unusual whirling arc in which a 
million volts jump between three points with vivid Hashes 
also leaping high in the air. 


* 
New Substage Lamp 


The Spencer Lens Company has recently developed a 
new microscope substage lamp which has a Bakelite 
molded housing. This housing is so strong that it may be 
knocked off the table and stood upon without breaking, 





Kinetic Theory of Gases. EarLeE H. KENNARD, professor 
of physics, Cornell University. First Edition. The McGraw- 


Hill Book Company, 
figs. Price, $5.00. 

The author remarks 
kinetic 


New York, 1938. xiii+483 pp. 94 
his preface that although the 
gases is a branch of physics which has 
passed from the stage of excitement and novelty into 
staid maturity, 


theory of 


there remains a demand for adequate 
exposition ef the subject. To meet this demand he has 
provided a clear and concise presentation of kinetic theory, 
from which the illogical and unsystematic elements have 
been excised so far as possible. While primarily classical 
in point of view, the book treats in detail certain reinter- 
kinetic theory quantum 
These particular 
significance in connection with collision processes, specific 


yretations of demanded by 
I . 


mechanics. reinterpretations are of 
heats, and the properties of the ideal gas from the stand- 
point of the Bose-Einstein and Fermi Dirac statistics. 
After an introductory chapter devoted to 
and fundamental concepts, 
are derived 


definitions 
the general distribution laws 
and discussed by conventional methods. In a 
chapter on the spatial distribution of molecules there is a 
thorough treatment of 
wave mechanical 


from the 
as well as from the classical point of 
The middle portion of the book is chiefly concerned 
with the theory of transport 


molecular scattering 
view. 
processes, the equation of 
state, and the thermodynamic functions of real and ideal 
gases. An interesting chapter on statistical mechanics, 
in which its foundations are critically discussed, is followed 
by a detailed account of the wave mechanics of gases. 
The book concludes with a chapter on the electric and 
magnetic properties of gases. 

The author is to be commended for his competent and 


thorough treatment of a field, which, though not new, 
still presents many interesting problems. The book should 


prove particularly useful as a text in kinetic theory for 

seniors and graduate students engaged in the study of 

physics and chemistry. Joun G. KirKwoop 
University of Chicago 
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Eisenlose Drosselspulen (Air Core Coils). J. 
Pp. 316+viii. Figs. 253, Tables 32, 163233 cm. 
Koehler Verlag, Leipzig, 1937. 


It has long been a source of annoyance in the design of 


Hak° 
K. F, 


air core coils for various purposes to find the pertaining 
information widely scattered over many periodicals. The 
above book not only represents a very intelligent collection 
of practically all the formulas developed since Maxwell's 
famous treatise on Flectricity and Magnetism, but it also 
gives a critical survey of their utility emphasized by many 
numerical tables for the useful ranges. All the formulas are 
carried back to their origin and written in a consistent 
uniform style, are illustrated by 65 examples chosen from 
all fields of application and presented frequently in graph 
form. As a novel feature, all the figures are collected at the 
end of the text and furnished with English captions. In 
addition, a German-English glossary of the more important 
terms adds to the usefulness for readers not quite familiar 
with the German technical expressions. 

The arrangement of the text is logical, starting with the 
computation of self-inductance of circular and polygonal 
coils based upon the Neumann formula. The influence of 
insulation is briefly covered, as well as the important aspect 
of the most economical design. This is followed by the more 
complex computation of mutual inductance bet weens coils 
of circular, rectangular, polygonal cross sections and vari- 
ous relative positions. Mechanical forces, effect of heating, 


skin effect and distributed capacity are well taken into 


account in separate chapters and many special coils are 


investigated, as for example, the spherical coil, coils of 


stranded wire, ring coils, and others; examples of con- 


. Finally, 


are briefly summarized and in an appendix the most recent 


struction are illustrated by pictures tests on coils 
form of the dust core coils is described. 
best 
bibliography containing 687 references covering the litera- 


The thoroughness 


of this monograph is illustrated by the extensive 


ture up to early 1937. The book should prove most useful, 
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‘ndeed indispensable, for all those who have to design and 
use coils for common as well as special purposes. 

ERNST WEBER 

Polytechnic Institute of Brooklyn 


Booklets Recently Received 


The Ohio Chemical Manufacturing Company of Cleve- 
land, Ohio, has issued a 16 page booklet dealing with 
equipment for control and storage of high pressure gases. 

Metallographic Equipment and Accessories, Catalog 
E-225, published by Bausch and Lomb Optical Company 


of Rochester, New York. This 32 page booklet contains 
descriptions of Ils Metallographic Equipment, Si Metallo- 
graphic Equipment, and Vertical I]luminator. 

The General Radio Experimenter. (8 pp.) Vol. XII, 
No. 10, March, 1938. This issue contains articles on 
“Modern Improvements for the Beat-Frequency Oscil- 
lator,”” ‘‘Variacs Used as Controls in Flexible Stage- 
Lighting Unit,”’ and ‘“‘Measurng 100,000 Rpm.”’ 

Oscillographer. (8 pp.) Vol. 2, Nos. 1-2, March-April, 
1938. In this issue are articles on ‘Type 169 Nine 
Inch Laboratory Oscillograph,” and ‘The DuMont 
Simplified Television System Demonstrated.”’ 











Dr. William T. Richards 
received his A.B. degree 
from Harvard in 1921 and 
his Ph.D. from there in 
1924. He was assistant to 
the director of Gibbs Me- 
morial Laboratory from 
1923 to 1924. He held an 
International Education 
Board fellowship from 1924 


to 1926. He was an instruc- 





tor at Princeton from 1926 
to 1928 and later was ap- 
pointed assistant professor. Dr. Richards is now with the 
Rockefeller Institute for Medical Research in New York. 


Dr. Ernest G. Linder received his B.A. degree in 
1925, and his M.S. in 1927 from the State University 
of Iowa. In 1931 he received 
his Ph.D. from Cornell Uni- 
versity. He was an instructor 
in physics at the University of 
lowa from 1925 to 1927, and at 
California Institute of Tech- 
nology from 1927 to 1928. From 
1928 to 1932 he was a Detroit 
Edison Research Associate at 
Cornell. Since 1932 he has had 
the position of research physicist, 
Research Division, RCA Manu- 
facturing Co., Inc. 





Dr. Hans Klemperer received his doctor's degree in 
electrical engineering from the Institute of Technology in 
Berlin in 1929, He was an instructor at the University of 
Aachen until 1930 when he accepted a position in the 
research laboratories of Siemens and Halske in Berlin. In 
1933 he worked in the development department of the 
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Siemens Schuckert Werke in Berlin, in charge of applica- 
tions of grid controlled rectifiers. In 1934 he was connected 
with the Radio Corporation in the electronic tube labora- 
tory at Harrison, New Jersey. Since 1935 he has been with 
the Research Laboratory of Westinghouse Electric and 
Manufacturing Company at East Pittsburgh, Pennsyl- 
vania, working on ignitron and gaseous tube development. 


Mr. R. Weller received a 
B.S. from M.I.T. in 1927 
and subsequently taught 
mechanical subjects at Me- 
chanics Institute in Roches- 
ter, New York. He obtained 
an M.S. in physics from 
Ohio State University in 
1934 and has been engaged 
in photoelastic research in 
the department of physics 
and in the Engineering Ex- 
periment Station of Ohio 
State for several years. 





Dr. George H. Shortley received his bachelor’s de- 
gree in electrical engineering at the University of Minne- 
sota in 1930 and his Ph.D. 
in physics from Princeton in 
1933. He was subsequently 
a National Research Fellow 
at M.I.T. and Harvard and 
is now assistant professor of 
physics at Ohio State 
University. His principal re- 
search has been in the field of 
the theory of atomic spectra 
and he is co-author of the book 
Theory of Atomic Spectra, by 
Condon and Shortley. 
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Contributed Original Research 





Contribution to the Problem of Arcback in High Voltage Discharge Devices 


H. KLEMPERER 
Westinghouse Electric and Manufacturing Company, East Pittsburgh, Pennsylvania 
(Received February 4, 1938) 


It has been found in rectifier and ignitron practice that 
arcback can occur after deionization of the discharge path 
at such vapor pressures where no self-sustaining glow can 
develop. Experiments to study this kind of an arc break- 
down were made in Hg vapor vessels between cold elec- 
trodes at voltages between 10 and 60 kv. It appeared that 
the start of arc breakdown was not dependent on the 
Paschen law. At the voltages used and distances between 
the electrodes of less than 50 cm arc breakdown always 
started at PD values which were much lower than those 
required to sustain a glow discharge. This is contrary to 
earlier experience at lower voltages. With rising voltage 


DESCRIPTION OF PROBLEM AND EARLIER 
EXPERIENCE 


HE phenomenon which is commonly called 

“arcback” and which will be called ‘‘arc 
breakdown” in some of the following discussions, 
is noticed as a sudden appearance of an arc spot 
on the ‘‘anode’"! of a rectifier during the inverse 
part of the cycle, i.e. when the ‘‘anode”’ is nega- 
tive. The discharge which is observed following 
breakdown shows the typical properties of the 
arc discharge: low arc drop and high current 
density in front of the cathode. 

No direct experiments showing onset of arc 
breakdown have been made so far. Cathode-ray 
oscillograms* which, due to inertia of the trigger- 
ing relay (Norinder type) started after elapse of 3 
microseconds, did not reveal any sign of a pre- 
ceding phase of the discharge, during which 


Anode” and ‘‘cathode” in quotation marks mean 
mechanical parts following rectifier practice, and do not 
denote the electrical polarity. 

? Oscillographic investigation by Maxfield, published by 
Slepian and Ludwig, Trans. A. I. E. E. 51, 92 (1932). For 
oscillograph see O. Ackermann, Trans. A. I. E. E. 49, 467 

1930). 
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arc breakdown frequency increased in an exponential 
manner. Changes in vapor pressure and electrode spacing 
had no definite effect on arc breakdown frequency as long 
as glow discharge conditions were not approached. Changes 
in partial air pressure did not affect arc breakdown fre. 
quency if air pressure was below 10 percent of mercury 
vapor pressure. Heat treatment of the cathode (graphite 
‘‘anode’’) has a large effect in reducing arc breakdown. 
Even well heat-treated ‘‘anodes’’ show seasoning effect, i.e, 
steady decrease of arc breakdown frequency after initia] 
maximum. Phenomena on cathode surface are believed to 
be responsible for formation of arc breakdown. 


voltage drop was higher.* The highest speed 
type of the cathode-ray oscillograph has not yet 
been used for such investigation, as far as the 
author knows, and it may be expected that faster 
oscillograms would provide insight into the for- 
mation period of the arcback discharge. 
Though the problem of arcback is as old as 
the mercury arc rectifier itself, no complete 
theory has been formed so far and rectifier 
and tube designing to a large degree is still 
empirical. In early days poor vacuum, which 
caused the development of a glow discharge‘ 
during the inverse part of the cycle, was held 
responsible for the formation of arcback. The 
vacuum was improved with the development of 
better pumps and vacuum technique, and the 
number of arcbacks decreased. But they did not 
disappear. Later on ion diffusion to the ‘‘anode” 
was made responsible for arcback® and shielding 
of the anodes and deionizing grids greatly im- 
proved the quality of the rectifier. But backfires 


8 See Rogowski, Zeits. f. Physik 100, 1 (1936). 

4R. Seeliger, Naturwiss. 16, 655 (1928). 

> J. v. Issendorff, M. Schenkel and R. Seeliger, Wiss. 
Ver6ff. a.d. Siemens Werken, 9, 73 (1930). 
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were not eliminated. Then the ignitron was 
believed to be the ideal solution of the backfire 
problem," since it does not use ‘keep alive’ 
means during the inverse part of the cycle and 
ion diffusion from other arcs is eliminated. But 
even with the ignitron arcbacks were observed, 
and these arcbacks occur not only during the 
transition period, but also near the time of 
the inverse voltage peak, when residual ioniza- 
tion practically has disappeared. 

As mentioned, some arcbacks—mostly those 
in high current rectifiers—happen during the 
transition periods.’ These arcbacks apparently 
are associated with residual ionization of the 
discharge path, which was conducting a short 
time previously. Other arcbacks, particularly 
those encountered in high voltage rectifier and tube 
work, occur near the time of the inverse voltage 
peak with very little or no ionization present. 
It is difficult to connect this type of arcback to 
any such phenomenon’ as known to start a self- 
sustaining discharge. lon bombardment of the 
cathode is very small, thermionics are hardly 
involved in the initial phase, and the field 
strength at the cathode surface is of the order of 
1000 to 10,000 v. cm, which is very small com- 
pared to gradients at which cold emission is to 
be expected. 


RANGE OF TESTING, SET-UP AND PROCEDURE 


Breakdown between cold electrodes, with no 
voluntary source of ionization present, and more- 
over under voltage, vapor pressure and spacing 
conditions where neither the development of a 
(Townsend) glow nor cold emission is expected 
may have some other relation to field strength. 
With this in mind the author carried out a 
number of breakdown tests, applying high uni- 
directional voltage to electrodes of test vessels. 
These vessels contained a graphite ‘‘anode,”’ 
a mercury pool “‘cathode’’ and no voluntary 





®W. G. Dow, Fundamentals of Engineering Electronics 
(J. Wiley 1937), p. 488. J. Slepian and L. R. Ludwig, “A 
New Method for Initiating the Cathode of an Arc,” 
Trans. A. I. E. E. 51, 693 (1933). 

‘Ludwig, Maxfield and Toepfer, Trans. A. I. E. E. 
53, 75 (1934). 
_ *Sometimes radiation may have some influence; there 
is NO reason against the consideration of cosmic rays as a 
cause of some of the breakdowns. For other mechanism 
see reference 17. 
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source of ionization. Voltages between 15 and 60 
kv were used. The Hg vapor pressure was varied 
between 0.26u (5°C) and 89u (80°C) and electrode 
spacings ranged from 1 cm to 50 cm (sometimes 
100 cm). 

The circuit diagram of the test set-up is given 
in Fig. 1. Engineering experience seems to indi- 





+ ll0V 





Fic. 1. High voltage set. (0.2 should read 0.2uf; kv should 
be meter.) 


cate that high voltage systems with a limited 
power supply are less liable to show breakdowns 
than high power systems of same voltage and 
electrode spacing.’ The probable reason is that a 
low power system does not supply current 
sufficiently rapidly for the arc spot to develop 
after the breakdown is initiated, and that with 
some unnoticed glow or flash the failure may 
clear itself. Similar reasoning applies to low 
pressure breakdown, if the breakdown is detected 
by the following arc discharge current.” 

The current density in the cathode spot of an 
arc discharge lies around 5000 amperes, cm’. 
The cathode spot of the arc, formed by break- 
down as seen from visible marks on the ‘‘anode,”’ 
sometimes is as small as } mm in diameter. 
Therefore, the high voltage test set should be 
able to supply a current of 5 to 10 amperes 
during the first few microseconds. As shown by 
sketch the test set used a condenser discharge to 
feed the breakdown. Time constant was 7=50 
microsecond; initial current at lowest voltage 
(15 kv) was 30 amperes. 

A sketch of the test vessels used is given in 
Fig. 2. Vessel A and B were water-cooled steel 
tanks with fixed electrode distances. Vessel C 


* Compare spacing in a 100 kv x-ray set with spacing 
used for power stations of same voltage. 

1° In high voltage rectifier and tube work the number of 
arcbacks observed depends on sensitivity of back current 
measuring device. 
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Fic. 2. Test vessels. Dimensions in cm. 


had glass walls and the electrode spacing was 
variable between 1 and 50 cm. 

Testing was done in the following way: the 
vessel was evacuated to a certain residual gas 
which read with the McLeod 
gauge. The temperature of the vessel was held 
constant at a desired point by regulating the 
cooling water or using a furnace. When condi- 
tions seemed to be constant direct voltage was 
applied in the direction prevailing during arc- 
back periods (‘‘anode”’ negative) and the occur- 
rence of breakdown was watched or counted over 
a certain length of time. Records of the break- 
down number 


pressure was 


an electronic 
counter. These tests were repeated for a great 
number of times in the hope of obtaining repro- 
ducible average results. 


were made by 


DISCUSSION OF RESULTS 
(a) PD law not applying to breakdown 


Right from the beginning of the tests it 
appeared that the pressure distance relation 
(Paschen law or PD law), which states that a 
discharge starts at a constant voltage if the 
product of gas pressure Xelectrode spacing (PD) 
is kept constant, does not apply to arc break- 
down. With fair accuracy a glow discharge could 
repeatedly be started at the same voltage if 
pressure and distance were changed as the law 
prescribes. The start of arc breakdown, however, 
did not follow any such relation. 

Figure 3 represents this situation in the voltage 


328 


voltage ignitrons, were arc. 
back was observed frequently 
near the peak of inverse 
voltage, and very seldom during the transition 
period. On the other hand, at lower voltage, 
smaller spacing (D=4 cm) and accordingly 
higher vapor pressures Marshall and Toepfer! 
observed glow discharges before any breakdown 
started. Some of their results are given in the 
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lower part of the VPD diagram of Fig. 3. 
This situation seems to apply to some arcback 
phenomena in regular power rectifiers, which 
run at large current and moderate voltage 
(600 v for instance). Here, in many cases, arcback 
was notably preceded by a glow discharge. 


1 [t is more exact to say the ‘‘u=1 curve.” u is the rise in 
ionization, defined. u=y(e*“—1). See Rogowski, Zeits. f. 
Physik 100, 15 (1936). 

12 Smede, Electric J. 25, 403 (1928). 

13 Marshall and Toepfer, test results, partially pub- 
lished by Slepian and Ludwig, Elect. Engineering 50, 793 
(1931). 
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(b) Influence of voltage on breakdown frequency 


The well-known increase in breakdown fre- 
quency with rising voltage appeared from all 
experiments without exception. It was not 
possible to find a general mathematical relation 
because in most of the tests the probability of 
arc breakdown changed during the tests, even 
after hours of waiting. Especially if the voltage 
was high it was difficult to get repeatable average 
values. In general, the rise of backfire frequency 
with voltage had an exponential character,'4 
and in one case where a commercial high voltage 
rectifier was subjected to the described direct 
voltage tests, results checked fairly well with the 
law, as given by Fig. 4. In this test all ‘“‘anodes”’ 
of the rectifier were tied together to test voltage ; 
all grids were kept floating and “cathode” and 
tank were tied to ground. The rectifier was 
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Fic. 4. Anode voltage vs. breakdown frequency of 
commercial mercury arc rectifier. 


thoroughly treated and aged, and the above 
readings were taken after several hours of 
testing. 


(c) Influence of vapor pressure on breakdown 
frequency 


No definite tendency of the breakdown fre- 
quency to increase as a result of rising vapor 
pressure could be observed. If conditions were 
approached, under which a self-supporting glow 
discharge can exist, the breakdown frequency 
rose rapidly. But in the region below this happen- 
ing the backfire frequency drifted erratically up 
and down, and other uncontrolled influences 
overshadowed any effect of vapor pressure.'® 


4 Slepian and Ludwig, Trans. A. I. E. E. 51, 92 (1932): 
® This experience should not be applied to general 
rectifier practice without considering other effects of 
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Fig. 5 gives results of such measurements taken 
on vessel A. Similar observations were made with 
all the other test vessels. 
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Fic. 5. Breakdown frequency vs. vapor pressure at 60 kv 
d=22 cm (taken with vessel ‘‘A”’ on Nov. 27, 1935). 


(d) Influence of electrode spacing on breakdown 
frequency 


Similar erratic relations were observed be- 
tween electrode spacing and the breakdown 
frequency. If the geometrical configuration leaves 
unobstructed a path, along which a glow dis- 
charge can start, the backfire frequency in- 
creases. If such conditions exist, shields, baffles 
or grids which suppress the start of a glow from 
the upper side of the ‘‘anode” or from an un- 
shielded “‘anode”’ stem also reduce the frequency 
of breakdown. On the other hand, if the electrode 
spacing becomes small enough to raise the field 
strength up to cold emission limits, breakdown 
follows. This effect has been observed with rough, 
untreated ‘‘anodes’’ at distances around 1 cm 
and voltages as low as 50 kv. But between these 
two borders, glow and cold emission, no definite 
tendency of the breakdown frequency in relation 
to spacing could be observed. Fig. 6 gives a 
sketch of test vessel C. A metallic shield was 
placed inside the upper part of the glass cylinder 
and tied to ‘‘anode,’’ stem and cover plate. 
Lines of force are given on one side, breakdown 


increased vapor pressure: the increase in deionization 
time, the change of arc drop, and maybe increase in 
leakage rate. 
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Fic. 6. Breakdown frequency and lines of force in vessel 
C with metallic shield. Arrows indicate spots at which 
backfires occurred. Tests made at 50 kv, 2yu air, 30°C Hg. 


spots, as observed by experiment, are indicated 
on the other side. Apparently some of the break- 
down started from points with higher field 
strength, (edges), others developed along the 
line of the longest path. Repeated tests gave the 
same general situation but, of course, not exactly 
the same distribution of breakdown spots. 


(e) Influence of the partial air pressure on break- 
down frequency 


A great number of tests was run to find out 
the influence of air admixtures to mercury vapor 
on the breakdown frequency. Air was made 
leaking in from the atmosphere with a needle 
valve, and in other cases the slow development 
of a leakage pressure was waited for and meas- 
ured with the McLeod gauge. The results were 
equally indefinite and hard to repeat. If the per- 
centage of foreign gas was less than 10 percent 









of the mercury pressure no effect on breakdown 
frequency was observed. If the foreign gas 
pressure was higher than 10 percent the break- 
down frequency was sensitive to any change of 
the foreign gas pressure. Mostly the breakdown 
became more frequent with increase in admix- 
ture, but there were cases observed when the 
effect was reversed. The minimum breakdown 
voltage of air is about one-half of that of mercury 
vapor, and the corresponding PD value (330, 
cm) is ten times lower than that of Hg, so that, 
as far as glow is concerned, we should expect a 
considerable influence of small air admixtures. 
But for the breakdown frequency such 
influence of partial air’ pressure can be stated 
from the described tests. 


no 


(f) Influence of heat treatment and aging on 
breakdown frequency 


It is known from practice that a well-designed 
rectifier with good vacuum inside would show 
poor arcback qualities if untreated anodes were 
used. On the other hand, poorly designed recti- 
fiers have been brought up to a fair performance 
simply by subjecting the anodes to a severe heat 
treatment and aging process. 

Heat treatment consists in running the anodes 
up to 1400 or 1800°C in vacuum. Even after 
such treatments new anodes show initial back- 
fires which become more and more seldom in 
frequency with the time of use, provided that 
no damage is done to anode surface or vacuum. 
Untreated anodes show this ‘‘seasoning”’ effect to 
a much larger degree and do not reach the final 
quality of heat treated anodes. 

Both procedures, heat treatment and aging, 
evidently do the same thing: they remove back- 


fire causes from the ‘‘anode”’ surfaces. Whether 
these causes consist of thin local dirt layers,” 
gas occlusions, irregularities in surface structure 
or loose particles, is still a matter of speculation. 


16 For the effects of small H» percentages see Dallenbach 
and Gerecke, E. T. Z. 57, 937 (1936). ; 

7 A. Guentherschultze and H. Schnitger, Zeits. f. 
Physik, 102, 34 (1936). Malter, Phys. Rev. 50, 48 (1936). 
The possible relation of this type of low voltage-cold 
emission discharges across thin layers (‘‘jet discharges”) 
to the formation of arc breakdown will be discussed later. 
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Effect of High Energy Electron Random Motion Upon the Shape of the 
Magnetron Cut-Off Curve 


E. G. LINDER 
RCA Manufacturing Co., Inc., Camden, N. J. 
(Received March 2, 1938) 


The shape of the cut-off curve is explained on the basis of a component of high energy random 
motion superimposed on the electron orbital motion (see references). Electron temperatures 
determined from the cut-off curve shape agree with those determined by probe methods. The 
possible effects of the random motion on other magnetron characteristics are pointed out. 


GREAT many experimental investigations 

have been made of the shape of the mag- 
netron* cut-off curve, but no one appears to 
have attempted to determine what its shape 
should be from theoretical considerations. How- 
ever, Hull's! work on the motion of electrons in 
magnetrons, and some experimental data on the 
effects of space charge, permit certain qualitative 
conclusions to be drawn immediately. 

In a perfectly symmetrical cylindrical mag- 
netron, neglecting space charge effects and the 
initial velocity of the electrons, and assuming 
that the electrons move only under the action 
of the electric and magnetic fields as described 
by Hull, the cut-off curve would have the general 
shape shown in Fig. 1(a). The anode current J, 
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Fic. 1. Idealized magnetron cut-off curves. (a) With no 
space charge or initial velocity; (b) Space charge but no 
initial velocity; (c) With both space charge and initial 
velocity. 


remains constant until the magnetic field reaches 
the critical value J7, when J, drops abruptly to 
zero. This is evident from the equations of the 
electron path. If space charge is considered, J,, no 
longer remains constant as JZ] becomes greater. 
The increasing curvature of the electron tra- 
jectories, and the longer transit time, cause an 
increase in space charge and a consequent de- 


* The type of magnetron considered in this paper con- 
sists of a cylindrical anode having a concentric filamentary 
cathode, aligned axially in a uniform magnetic field. 

1A. W. Hull, Phys. Rev. 18, 31 (1921). 
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crease in J, before the cut-off point is reached. 
However, at J/=I/1/,, I, drops abruptly to zero 
as before, as is shown in Fig. 1(b). If initial 
velocities also are considered the curve takes 
the form shown in Fig. 1(c). The first part of 
the curve, for J7<JI,, is the same as in the 
previous case. But at /7=//, an abrupt drop no 
longer occurs. At that point only those electrons 
having zero initial velocity fail to reach. the 
anode, those having higher initial velocities 
require fields greater than /7, for cut-off. This 
yields a curve of the form shown, which is the 
form always observed experimentally if the 
emission current is sufficiently large to cause 
appreciable space charge. 

The first part of the cut-off curve, i.e., /7 <JI,, 
has been discussed previously by L. Tonks,? and 
more recently in greater detail by G. Griinberg 
and V. Wolkenstein,? and will not be dealt with 
in the present article which will discuss only the 
second part, i.e., J/7>JI,., the shape of which 
depends upon initial velocities. 

Previous work‘: ®> has shown that electrons in 
magnetrons may acquire a Maxwell-Boltzmann 
distribution of velocities superimposed upon their 
orbital motion, and theretore behave very much 
as though emitted from a cathode of very high 
temperature. Probe measurements have shown 
that these temperatures are many times that of 
the actual cathode being of the order of 40,000°C 
for a plate potential of 100 volts in tubes in- 
vestigated by the writer. Initial velocities corre- 
sponding to such temperatures should affect very 


2 L. Tonks, Physik. Zeits. Sowjetunion 8, 572 (1936). 

*G, Griinberg and V. Wolkenstein, Tech. Phys. USSR 
4, 179 (1937). 

‘IT. Langmuir, Phys. Rev. 26, 585 (1925). 

*E. G. Linder, Proc. I. R. E. 26, 346 (1938). 
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considerably the electron orbit and the shape of 
the cut-off curve. 

An expression will be derived giving J, as a 
function of JJ and the temperature 7° of the 
electron swarm in the cathode neighborhood. 
In order to derive this expression, assume that 
all electrons start so near the central axis of 
the tube that only radial and axial components 
of initial velocity need be considered, components 
tangential to the cathode being negligible. As- 
sume also that the initial velocities correspond 
to the temperature 7 of the Maxwellian dis- 
tribution. 

The plate current due to electrons of number 
dn having initial radial velocities between u and 
u+du is given by 


dI,=ev dn, (1) 


where v is the radial velocity at the plate. 
From Hull,! the radial velocity at any distance 
r from the origin is given by 


2eV Tle\? i 
- r+ul, 
m 2m 
where V is the potential at any point, and J/ is 
the magnetic field strength. 
The number of electrons having initial radial 


velocity components between u and u+du ac- 
cording to Maxwell's distribution law, is given by 


dn = 2nye~**x dx, (3) 


where x=(m 2kT)*u, and ny is the total number 
emitted. The need not be 
considered since they do not contribute to the 
plate current. 

The cut-off value J/, of the magnetic field is 
found by putting dr/dt=0, for r=r,, and u=0, 
in (2), 


axial components 


i.e. JI,.=(8mV’,/er,7)?. (4) 
Introducing (2), (3) and (4) in (1), yields 
dI,=A(x*—a*)'e x dx, 
where a?=(re/2m)*m 2kT (UI? —H,,). 


We are concerned only with //>//,, i.e., with 
positive values of a*. Hence J, will be real only 
for x*>a*. Hence we integrate from a to «, 
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D 


1,=A | (x? —a*)*e-**x dx, 


a 
T,=Ts a? — Ji_(m/2kT)(re 2m)?(H ¢? H*) | 


9 


I, m fre? 
log = ( ) (11,2 — HH). (5) 
I, 2kT\2m 


J, in this expression is considered constant, 
its variation compared to that 
of J,. J, is the total current emitted by the 
cathode—not the real cathode but the high 
temperature virtual cathode discussed previously 
by the writer.® At the cut-off point, i.e., J7=H,, 
it is equal to the space charge limited current J,, 
and may be computed from the Langmuir-Child 
law, proper correction being made for electron 
trajectory curvature due to the magnetic field. 
For JJ>JI,, the returning current due to those 
electrons which do not reach the plate and 
return to the cathode must be considered. This 


being small 


returning current J. plus the advancing current 
7; must approximately equal the space charge 
limited current, i.e., J:+/.=Jo. For complete 
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Fic, 2, Log J, versus H? from Hull’s data (reference | 
Fig. 8, page 44. 
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Fic. 3. Log J, versus H® showing linear portion of curve 
and change of slope with plate voltage. Curves are labeled 
with plate voltage V>. 


cut-off the plate current becomes zero and we 
have the returning current equal to the ad- 
vancing current, i.e., Jj;=J.=J)/2. It is thus 
seen that Jp/2<1,< Io, and there is a maximum 
variation of J;, with respect to J» of a factor 
of two. 

The space charge limited current J itself also 
may vary as // is increased, due to changes in 
the size of the electron swarm, and potential 
distribution. The magnitude of this change may 
be estimated by taking the Langmuir-Child 


equation 
2/2e\} Vi} 
9\m r 


where V; is the potential at the edge of the swarm 
and r; the radius of the swarm. The magnetic 
field may be expressed in terms of V; and r; by 
an expression similar to (4) 


i.€., IT=(8mV,/er,?)}. 
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Eliminating r,; between these two equations 
yields Iy>=(e, 9m) V\1I, which gives the variation 
of the space charge limited current with J/ 
and JV. As JI is increased the swarm shrinks in 
size and therefore V, decreases. Thus the 
changes in JJ and VV, are compensatory. But 
even without the compensating change in |’, the 
change in J») would be negligible compared to 
that of J,. 

The total change in J, is thus seen to be small 
compared to the change in J,, and J, will there- 
fore be assumed constant in discussing the 
experimental results. 

From (5) it is evident, that if log J, be plotted 
against J/* a straight line results whose slope 
depends only upon 7. Data from Hull’s paper 
plotted in such a fashion are shown in Fig. 2. 
Similar data obtained by the writer are plotted 
in Fig. 3. The portions of the curves to which 
the above analysis applies are fairly linear. 

It is evident that the slope, i.e., electron tem- 
perature, increases with the plate potential V,. 
Data from Figs. 2 and 3 showing this are plotted 
in Fig. 4. The points are somewhat scattered, 
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Fic. 4. Electron swarm temperature in equivalent elec- 
tron volts versus plate voltage, as determined from cut-off 
curve and by probe measurements. 


but within the experimental error lie on a straight 
line through the origin. This is in agreement 
with previous measurements of electron tempera- 
tures by probe methods,*: ° in which the electron 
temperature was found to be proportional to V’, 
under similar conditions. 

Probe measurements were made in the present 
case also. These revealed that the electron tem- 
perature did not remain constant as // varied, 
but increased as cut-off became more complete. 
The range of variation is shown by the vertical 
arrow in Fig. 4. It may be seen that the mean 
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value given by the probe measurements is in 
fair agreement with the values determined from 
the cut-off curves. That the latter points seem 
somewhat too high may be due to tube dis- 
symmetries which would cause the cut-off curve 
to have a greater slope, and cause the tempera- 
ture determined therefrom to be too high. 

It may be concluded that the shape of the 





cut-off curve can be satisfactorily explained on 
the basis of the random electron motion. The 


effect of this motion on other magnetron charac- 
teristics such as orbit shape, transit time, transit 
angle, tube noise, etc., appears not to have 
been considered in any publications up to the 
present time, but there can be little doubt that 
these are all affected in an important manner. 





The Numerical Solution of Laplace’s Equation 


G. H. SHORTLEY AND R. WELLER 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
(Received February 28, 1938) 


This paper considers in detail numerical methods of 
solving Laplace’s equation in an arbitrary two-dimensional 
region with given boundary values. The methods involve 
the solution of approximating difference equations by 
iterative procedures. Modifications of the standard Lieb- 
mann procedure are developed which lead to a great in- 
crease in the convenience and rapidity of obtaining such a 
numerical solution. These modifications involve the use of 
formulas which simultaneously improve a block of points 
in place of a single point; methods of operating on the dif- 


ferences of trial functions in place of the functions them- 
selves; and also a method of extrapolating to the final 
solution of the difference equations. The theory underlying 
these procedures is considered in detail by a new method 
which involves the expansion of the error and difference 
functions in terms of eigenfunctions. This permits definite 
comparison of rates of convergence of various procedures, 
The techniques of handling practical problems are con- 
sidered in detail. 


$1. INTRODUCTION 


HERE are many industrial and engineering 

problems in which solutions of Laplace’s 
equation are needed. These solutions cannot 
often be obtained analytically because of the 
occurrence of irregular boundaries for which 
systems of separable coordinates are unknown. 
For example, our interest in this equation was 
aroused by the need for a convenient method of 
obtaining from the fringe photographs of photo- 
elasticity the principal stresses in the interior of 
the specimen tested. The photographs give 
directly only the difference of the two principal 
stresses. On the boundary of the specimen the 
normal component of the stress vanishes, so 
here we have complete information. The sum of 
the principal stresses satisfies Laplace’s equa- 
tion; values of this sum in the interior of the 
specimen can be obtained from the fringe photo- 


graphs by solving the two-dimensional Laplace 
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equation, using the known boundary values.' A 
numerical method must be used to obtain this 
solution since the shape of the specimen is 
usually much too irregular to admit of analytical 
treatment. The same necessity for numerical 
solution of Laplace’s equation arises in the cal- 
culation of electrostatic fields in regions enclosed 
by boundaries at known potentials; in the cal- 
culation of steady state temperature distribu- 
tions in heat conduction problems; in the cal- 
culation of the shapes of films and membranes; 
and in certain problems of hydrodynamics and 
gravitation. 

We shall restrict our discussion to the solution 
of the Laplace equation in a two-dimensional 

! Details of the application to photoelasticity of the 
theory outlined in this paper will be found in Weller and 
Shortley, ‘Calculation of Internal Stresses in Photo- 
elasticity,”’ to be published soon in the Journal of Applied 


Mechanics. This reference should be consulted for an 
explicitly solved example. 
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region with given arbitrary boundary values, 
with some mention of the case where the normal 
derivative at the boundary is given instead of 
the boundary value. The methods we shall give 
can be readily extended to three dimensions, but 
the amount of labor involved in solving a three- 
dimensional problem will be tremendously 
greater than that required for a similar two- 
dimensional problem. 

Any numerical method of handling a differ- 
ential equation necessarily involves the replace- 
ment of the continuum of points on the boundary 
and in the interior of the region by a discrete set 
of points. The best method which has been 
previously given for handling Laplace’s equation 
is an iterative procedure due to Liebmann.’ In 
this, a square net is laid down over the region; 
approximate values are assigned to the function 
at the interior points—known values at the 
boundary points. One then traverses the net 
repeatedly, replacing the value at each interior 
point by the sum of the values at the four 
neighboring points, using the new values im- 
mediately in the improvement of the succeeding 
points. The function will then converge to a 
solution of the Laplace difference equation 


(w(x, y+h) +(x, y—h) +0(x+h, y) 
+w(x—h, y) —4w(x, y)]/h?=0. (1) 


which will be an approximation to the solution 
of the Laplace differential equation—the smaller 
the interval h, the better the approximation. 


2 An examination of the literature of mathematical and 
applied physics indicates that numerical methods of 
solving Laplace’s equation have received little attention. 
Liebmann’s paper appears in Sitzungsber. der math.-phys. 
Klasse der Bayer. Akad., Miinchen, 1918, p. 385. The 
method is discussed by Phillips and Wiener, J. Math. and 
Phys. 2, 105 (1923); Courant, Zeits. f. Angew. Math. u. 
Mech. 6, 322 (1926); Frank von Mises, Differential und 
Integralgleichungen der Mathematischen Physik 1, p. 735. 
In these references little attention is given to the applica- 
tion of the method to practical problems: none of the 
material contained in the present paper is discussed. A 
different successive approximation method for the solution 
of the difference equation is given by L. F. Richardson, 
Trans. Roy. Soc. A210, 307 (1910); Math. Gazette 12, 
415 (1925). This method is similar to the Liebmann pro- 
cedure but does not use the new values, as they are ob- 
tained, in the improvement of succeeding points. This 
makes the convergence much slower (cf. section e., 
Appendix) unless considerable care is exercised in the 
assignment of certain arbitrary multiplying factors. The 
method does not have the general utility of the Liebmann, 
and at its best does not converge as rapidly as the modified 
methods we shall propose. 


VOLUME 9, MAY, 1938 





It has been shown quite generally by Richard- 
son* that the difference between the solution of 
a differential equation and that of an approxi- 
mating difference equation on a net of interval 
h is of the form 


A(x, y)h®+ B(x, y)ht+C(x, v)h® +--+. (2) 


Since only even powers of / enter here, the dif- 
ference solution approaches the differential 
solution rapidly as # approaches zero. If h is 
sufficiently small the error is proportional to h*, 
and by making two solutions with different h, 
one can estimate the difference of each from the 
solution of the continuous problem. 

In this paper we shall first investigate, by a 
new method, the rate of convergence of the 
iterative process to the final solution of the dif- 
ference equation in the case of the Liebmann 
procedure. These considerations lead to the 
formulation of procedures of similar type (in 
particular the “‘nine-block”’ procedure) which 
converge much more rapidly than the Liebmann. 
They also lead to the development of an accurate 
method for extrapolating to the final solution of 
the difference equations. In this way we greatly 
reduce the amount of labor required in obtaining 
a numerical solution of Laplace’s equation to any 
desired accuracy and render the task of cal- 
culating the solution entirely practical and con- 
venient. The last section of this paper is devoted 
to a detailed discussion of practical techniques. 


§2. THrory 


In this section we shall for simplicity assume 
that we are working with a square net of interval 
h in a region with a rectangular boundary. The 
slight modifications of the theory needed in case 
of a more general net or an irregular boundary 
will be discussed in §3 under the heading of 
Technique. We shall first examine the rate of 
convergence of the Liebmann iteration procedure. 

Denote the values of a function at the interior 
points of the lattice by ¥;, where i=1, 2, ---, n. 
Denote the values at the boundary points by yu, 
where a=n+1, n+2, ---, m. For example, for 
the lattice of Fig. 1, n=54, m=88. 

In the Liebmann procedure for solving the 
Laplace difference equation (1) with assigned 
boundary values ¥, one starts by assuming any 
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Fic. 1. A net with 54 interior points 
and 34 boundary points. 


values Y; whatsoever at the interior points. Then 
one works over the lattice, e.g. in the order in 
which the points are numbered in Fig. 1, replac- 
ing each y; by the mean of the four neighboring 
values to get an improved value y;. The pre- 
viously improved values are used in calculating 
this mean wherever they are available, so that e.g. 
in Fig. 1, 


¥is=4}(Vi2tWitvistye). (3) 


The function y is then closer to the solution of 
the difference equation than was y.* Iteration of 
this process leads to values which are as close to 
the true solution as may be desired. 

In this improvement,‘ ~; may be expressed as 
a linear combination of all the y; and the yz, 
since for example in (3) the p's on the right are 
themselves combinations of other y's and y's 
which in the last analysis reduce to combinations 
of y's alone. Hence we may write 


Vi an > Tiibj+ } 3 U iaba, 


Ya=Ya, (4) 
where the 7;; and U;, are numerical coefficients 
which depend on the size of the net and the 
order of improving the points, but not on the 
particular boundary values given nor the par- 
ticular body values assumed. 

Now let w represent the true solution of the 
difference equation (1). The ¢rial function y 


3 See the convergence proof in the Appendix, or in the 
references noted in footnote 2. 

‘The arguments of the next four paragraphs apply 
equally well to the other types of improvement procedures 
considered in the balance of the paper and in the Appendix. 
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may then be expressed as the sum of the frye 
solution w and an error ¢: 


Vi=ait+ di, Ya= wat ba= wa. (5) 


Since the trial function must be chosen with the 
correct boundary values wa, the error ¢ will have 
vanishing boundary values: ¢.=0. Since the 
transformation (4) is linear, the function ob- 
tained when y is “improved”’ will be the sum of 
the functions obtained when w and @¢ are im- 
proved: ~Y=a+¢. But G=w—the improvement 
process is such that it does not change the true 
solution—this is seen to be true when the im- 
provement formula (3) operates on any solution 


of (1). Hence 
v=wt9, (6) 


and in order to study the rate of convergence of an 
arbitrary function Y with correct boundary values 
to the true solution w on iteration of the improve- 
ment, we need only investigate the rate of convergence 
of an arbitrary “‘error’’ function with zero 
boundary values to zero as the same improvement 
formula is used repeatedly. The two rates of con- 
vergence are the same. Since ¢,.=0, we have 
simply 


oi= LT iid; (7) 


or briefly 
$=T9. (7') 


Now the possible functions ¢ are best classified 
in terms of the eigenfunctions (characteristic 
vectors) of the matrix 7. There will be a set of 
functions ¢°, ¢@, , @*” such that 


To =h,o™, 


where \; is a constant. We show in the Appendix 
that all of the |X\'s will be less than unity. This 
set of eigenfunctions may be considered as form- 
ing a complete set® in terms of which any function 


5 The arguments of the text at this point are not rigor- 
ously true and therefore the conclusions drawn from them 
need justification which is furnished in this footnote. The 
discussion in the text would be mathematically correct 
only if the matrix 7 were symmetric, which a calculation 
for a simple lattice shows to be not the case. For example 
for a square lattice of four interior points, improved in the 
order 

1 2 
3.4 
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@ (with zero boundary values) may be expanded : 
6=a19" +429 + ---+a,90™. (13) 


the matrix T is the following nonsymmetric array: 


0 8 8 O 
0228 
0228 (8) 
0114 


(Here the elements of the top row are Ti, Tix, Tie, Tis; 
etc.) A nonsymmetric matrix cannot in general be reduced 
to diagonal form; its roots will not necessarily be real, and 
it will not have a sufficient number of eigenvectors to permit 
expansion in terms of them. But it can be reduced to a 
unique stepmatrix with steps of the form 


» 10 - 0 0 
0» 1 - 0 0 
. . : . . ° (9) 
000. rx 1 
000-02 


Associated with the root A, there will hence be a group of 
invariant vectors x‘2), x‘), --+ satisfying the relations 
Tx(™ = gx (a) 
5 Natit aR x, ed (10) 
T x6 = dgx( + x(a) 
of these only xm is an eigenvector. The entire set of 
invariant vectors is complete and the expansion (13) must 
I 


be made in terms of them. The canonical form of the 
matrix (8) is 


1 0 0 0 
000 0 
0001 (11) 
0000 


while the canonical form of the matrix for a six-point net 
improved in the order 
i223 
e 3 
3643 0 0 0 O O 
0 0625 0 0 O 
0 O 01070 0 O 
0 0 0 0 9O 
0 0 0 0 =O 1 
0 O 0 0 OO O 


In spite, however, of the fact that T is not a symmetric 
matrix, it is convenient to speak in the text as if it were, 
because of the simplicity and familiarity of the simple 
eigenvalue theory of symmetric matrices. This is permis- 
sible because we are not led into any essential error thereby. 
It may be shown that repeated application of a stepmatrix 
composed of steps such as (9) to an arbitrary vector re- 
duces finally to the eigenvector of the root of highest 
absolute value. Thus the conclusions of the text are correct. 
The rapidity of this convergence to the eigenvector of the 
highest root is however less, the greater the size of the 
step involving this particular root. But it seems likely from 
geometrical and asymptotic considerations that the largest 
roots are all real and belong to steps of order 1, and that 
complex roots and invariant vectors which are not eigen- 
vectors occur only for the very small roots, where they 
are quickly wiped out by the iteration and do not cause 
trouble. These statements are borne out by. the form of 
(11) and (12), in which the only invariant vectors which 
are not eigenvectors belong to the root 0. Hence all the 
text arguments may be said to be ‘‘practically”’ true. 

For reference to the theory of nonsymmetric matrices 


see Bécher, Higher Algebra, Chap. XXI1, or Wedderburn, 
Matrices, Chap. IIl. 
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Then To=ia1¢6 +r2a2g + - - -, 
TTo= T26 =) 7019? +r22d2g + ve, 


Tb =)y'a1d™ +rs*aop™ + - - 


So after k applications of T, i.e., after k traverses, 
the coefficients of the various eigenfunctions in @ 
are decreased by \*. After a time only the function 
with the highest eigenvalue, say \,, will remain. 
The whole rate of convergence of ¢ to zero, and 
hence of y to w will depend on the rate at which 
the ¢“ part of ¢ is wiped out—i.e., on the size 
of 1. 

These considerations furnish a practical method 
of finding \; and ¢™ for a given lattice and order 
of improvement. ¢“ will clearly be “smooth,” 
nodeless, and will have a maximum near the 
center of the region in question. Take a trial 
function of this general shape, iterate repeatedly, 
and soon all the terms involving ¢’s of lesser 
eigenvalue will have dropped out, leaving a 
function proportional to ¢, which on each 
traverse is reduced uniformly to \, times itself.® 
In this way, the values of \, for regions with 
square boundaries of various sizes were deter- 
mined for the Liebmann procedure. The values 
are plotted in Fig. 2 in the form of a curve 
giving — p” log, A; as a function of p, the number 
of intervals on a side of the square [(p—1)?=n, 








6 The determination of this \; value by iteration is 
facilitated by the following theorems, w hich apply to the 
— ement of functions of zero boundary value. 
If the nowhere negative function goes on improvement 
Pe a function » such that for all lattice points w<¢i/di<v, 
then w<d:<v. This theorem follows if we show that the 
bounds on ¢/¢ imply the same bounds on T¢/¢, and this 
follows by noting that in (7) all the 7;;’s are nonnegative 
(see Appendix for discussion of this property of 7;; in con- 
nection with the various types of improvement formulae). 
2. When one has arrived at the stage where all the 
eigenfunctions have disappeared except the first and 
second, and if ¢ is everywhere positive, one will find that 
along the nodes of $, ¢;/¢;=A1 and is constant on itera- 
tion; in regions where ¢°) >0, ¢:/¢:< \; but increases and 
approaches A; with iteration; in regions where ¢ <0, 
¢i/¢i>1 but decreases and approaches \;. To prove these 
statements consider a function of the form ¢=¢+¢™. On 
iteration ¢~>o= T¢—T?¢. Investigation of the successive 
ratios, T¢/¢ and T¢/¢ gives the above results when we re- 
member that A2»<A:. Hence when one arrives at the point 
where the ratio ¢/ between successive functions is decreasing 
with iteration over roughly half the region, increasing over the 
other half, and standing constant along the line dividing these 
regions, the value of the ratio along this dividing line is \,. By 
making use of this theorem one need not wait for ¢ to 


disappear but only for ¢\), , to disappear in determin- 
ing Ai. 
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procedures. The numbers opposite the points are the values 
of \; for a square region with p intervals on a side. 


the total number of interior points |. This form 
was chosen because asymptotically for large p, 
— p* log \1—7”’, i.e. for large p 


—log \,=1—A,=7°/ p’. (14a)? 


7 In the case of a rectangular region of p by g intervals, 
this asymptotic formula becomes 


—log A= 1— A= Hat ( a :) ° (14b) 
| an 
A lower eigenvalue \,,,, which belongs to the function 
with u—1 nodes in one direction and v—1 nodes in the 
other, has the same asymptotic value as that for the node- 
less eigenfunction in a region of p/u by q/v intervals. 
Hence asymptotically 
2 p? 
—log Xy, >= 1—y, »= da? (5 + ») (14c) 
@? 

The derivation of these asymptotic relations is discussed 
in the Appendix. Although for small regions the form of 
¢” changes considerably with a change in the order of 
improvement of the points, the value of A; is to a good 
approximation independent of this order. The curve of 
Fig. 2 was obtained by the use of similar orders for the 
various size regions, but for all practical purposes it is 
valid for any order. No particular order seems to have ap- 
preciable advantage over any other. While we have no 
rigorous proof, we believe the asymptotic value to be 
independent of the order of improvement. 
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The number of traverses which must be made in 
order to reduce the coefficient of ¢™ in (13) to 
a definite fraction p of its initial value js 
(—log p)/(—log 1), so the value of —log d, is q 
measure of the rapidity of convergence of the process. 
Approximately, the rapidity of convergence 
varies inversely as the number of points in the 
square. 

We can then visualize the improvement process 
as follows: after a number of traverses have been 
made, the error will be a pillow-shaped multiple 
of @. A traverse consists in going over this 
pillow bringing successive points down to the 
mean level of their four neighbors. In this way 
the height of the whole pillow is gradually 
reduced to zero. We then see that if we could 
reduce at one step a greater area of the pillow 
to the ‘‘mean”’ level of its neighbors, the pillow 
would disappear faster—we would be cutting it 
down by larger slices. This more rapid slicing 
would be accomplished if we divided the interior 
of the region up into square blocks each contain- 
ing v points and then used improvement formulas 
which would make a whole block at a time a 
harmonic function [in the sense of (1)] with 
boundary values determined by the values of 
the function at the 4v points just outside the 
block. Such improvement formulas can be ob- 
tained by solving the v* simultaneous equations 
which state that the value at each point in the 
block equals the mean of its four neighbors. 

For example, we might divide the region of 
Fig. 1 up into blocks of nine points (nine-blocks) 
in the manner indicated in the drawing. We 
could then improve the six blocks successively 
in any order, for example the top row from left 
to right, then the bottom row from left to right. 
In this case, the improved values in the upper 
left-hand block would be calculated from the 
boundary values W356, Ws7, Wss, Ws, Wis, Wee, Wao, 
Woo, Vos, Ws6, ¥s7, Wss- Phe simplest formulas which 
give harmonic improved values are 
V2 = 41 2(37Ys2t+ 1 lvset Llysst+ 11 ysst+ lly, 

+72 t+ Tyist3ysot 322+ 3y2s+3ys0+5y29); 
vio, Viz, Veo are given by corresponding formulas; 
¥i=hVset¥sstvet+yio); (15) 


vs, Vig, ¥21 are given by corresponding formulas; 


vir = 1 (YotviotVizty¥20)- 
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In the improvement of the next block, 3, Pie, 
Yor would be used as three of the boundary 
values. 

By improving by blocks in this way one obtains 
considerably lower \ values and hence con- 
vergence to the final solution after fewer traverses 
than for the straight Liebmann process. We show 
in the Appendix that the asymptotic values of 
—p* log \ for a v*-block are v times the values 
given by (14), so the convergence is v times as 
fast as for the Liebmann. But the formulas used 
are so complex that it takes v times as much 
labor per traverse, so that little saving in work 
is made in this way. 


Nine-block procedure 


But it is possible to modify the improvement 
formulas* so that the improvement can be by 
nine-blocks with practically no more work per 
improvement than in the case of the Liebmann. 
In the nine-block 


x x x x 
ne ds e2 de ns 
x e 


x a ae ” 
e3 C3 ar C1 ea 
x 2 * x 
ds bs C4 b, dg 


x 
x 
x 
x 


we require that the value 


a= igLditds+d3+d,+d;,+d,+d;+ds) 
+2(e:+e2.+e3+e,) | 


as would be obtained by using (15); but we may 
now obtain the values b; by requiring that they 
be the means of their four diagonal neighbors, 
i.e. that 


* This modification amounts to a change in the set of 
difference equations which are satisfied by the final solu- 
tion, so that the final solution for this procedure will be 
slightly different from that for the Liebmann. That this 
does not in general impair the accuracy of approximation 
to the continuous solution is indicated by arguments in 
Section ¢ of the Appendix. 
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b,=4(ai+e:+e2+m), etc.,° 
and finally that the ¢ points be the means of 
their immediate neighbors: 


€,=}(a,;+),+e:4+5;),  etc., 


using in the last two formulas the improved a 
and b values. 

Precisely the same reasoning which we have 
applied to the Liebmann procedure may be used 
in the study of the rate of convergence of this 
nine-block procedure; empirical values of ¢ 
and X, for nets of various sizes may be similarly 
obtained, except that in this case the disap- 
pearance of all functions below ¢ will take 
place much more rapidly since all eigenvalues 
will be much smaller. The values of — p* log \, 
obtained in this way are plotted in Fig. 2—they 
are seen to be about 3.5 times the corresponding 
values for the Liebmann procedure. The asymp- 
totic value of —p* log \, for this procedure is 
ign" = 34.5 and other asymptotic eigenvalues for 
this procedure are given by writing 747? in 
place of 7 in (14) (see Appendix). These values 
are essentially independent of the order of im- 
proving the blocks. 

Since the number of traverses necessary to 
accomplish a certain reduction in the coefficient 
of an eigenfunction in the error is inversely pro- 
portional to —log \;, the number of times the 
Liebmann procedure must be iterated in order 
to accomplish reduction by the same factor as in 
one traverse by the nine-block procedure is just 
the ratio of the curves of Fig. 2. Thus we may 
conclude that in general the nine-block procedure 
will cut the time required to solve a problem to given 
accuracy by a factor of almost 3.5 as compared to 
the straight Liebmann method. This remains true 
when an extrapolation procedure to be dis- 
cussed later is used to remove the first eigen- 
function after the second has disappeared. The 
slightly greater inconvenience in improvement of 
the center-point of the nine-block is made up by 
a certain advantage in this extrapolation pro- 
cedure. 


* In case one of the » points is not available, the follow- 
ing formula, derived by the methods outlined in the 
Appendix, may be used for the corresponding 6 point: 


b, =} (2a, +e:+2d,+2d.+e2). 
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Four-block procedure 


It is possible to use blocks of four in place of blocks of 
nine, by the following formulas: 


x x x 
ta eo dg hy 
x € e x 
fea c2 a dy 
x « a x 
é2 6 Ci ex 
x x x x 
t3 gi fa ta 


a= log 7d + 7d:+2e,+2e:+2f;+2fo+2:+ 22), 
b=hlatfitfetis), 


C= i(at+b+e,+/f;), etc. 


The formulas should be used in this order, using the im- 
proved values of a and b as they are obtained. The value 
of \; varies according to which of the four points is taken 
as the point a to which to apply the long formula, in rela- 
tion to the order of improvement. The upper curve of Fig. 
2 is for the case in which the long formula is applied to the 
point of each four block which is nearest the boundary of 
the region, the lower curve when this formula is applied 
to the point nearest the center—improving the blocks row 
by row as one reads the words of a book. The two asymp- 
totic values given are for this same order of improvement, 
the higher one, 44; 97, when the long formula is applied 
to one of the points on the left side of each block, the 
lower, 27?, when this formula is applied to a point on the 
right side. 

One sees that the four-block converges somewhat more 
than twice as fast as the straight Liebmann, and is a 
definite saving in time, in spite of the fact that the labor 
per improvement is even more than that for the nine-block. 


Twenty-five-block procedure 


In the direction of blocks of size larger than nine, one 
could use a block of 16; in this case it would be necessary 
to find all four central values directly from the boundary 
values by means of fairly cumbersome formulae. One can 
do a twenty-five-block almost as easily, calculating only 
four points by cumbersome formulae (see Fig. 3). One 
first obtains the values at the points aj, de, a3, a4 by for- 
mulas similar to the following: 


a,>= lsog 38d, + 38d.+76¢:+ 76¢e>2 +49f,+49f.+26¢, +26, 
+12/,+12/.+6m,+6m.). 


Then the values at the points marked b are obtained from 
these and the boundary values by taking the mean of the 
four diagonal neighbors,’ and finally the values at the 
points marked c are obtained from the means of the four 
closest neighbors. This procedure gives an asymptotic 
value — p* log \, = 98447?=52.0, and hence considerably 
more rapid convergence than the nine-block. But it is 
probably quicker on the whole to introduce only nine- 
blocks since they are so much simpler and easier to improve. 
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One compensating factor in favor of the 25-block is that 
the central 6 and the four central c values need not be 
calculated until the iteration is completely finished since 
they do not enter into the improvement formula for any 
other points. , 

It is of course possible to handle still larger blocks by 
similar methods with even lower \,; values, but in these 
cases sO many points would need very complicated im- 
provement formulas that they would certainly not save 
time on the whole. The optimum block size seems to be 
nine, with perhaps some arguments in favor of 25. 


x x x x 
ne he £2 fe €2 da my 
e + e oe . 
t2 ba C3 be C1 bs dy 
7 + * * 7 ® x 
Je2 cs a2 C9 ai C2 ex 
e * a x 
ko by Cir bg Cio) Os bs fa 
* a ae a s x 
lo C7 a3 C12 aa C4 &. 
x a © * es x 
ma b.4 Ca bs Cé b3 hy 
2a x a x x x 
Na my ly ky ja ta N3 
Fic. 3. 


§3. TECHNIQUE 


In this section we shall summarize the prac- 
tical techniques which our experience has shown 
to be best adapted to the solution of the problem 
in hand. 


Choice of net 


The first thing of course is to lay over the 
region in question a lattice. It is usually most 
convenient to take a square lattice, although 
there are probably cases where a rectangular 
lattice is advisable, and cases where it is advisable 
to change the lattice spacing near sections of the 
boundary where the function representing the 
boundary values has exceptionally high or low 
curvature.” The net points at which the value 
of the function is to be defined are the interior 
points at the intersections of the bars of the 
lattice and the boundary points at the inter- 
section of the lattice bars with the boundary of 
the region, as indicated in Fig. 4. 


10See Appendix for improvement formulas for these 
cases. 


JOURNAL OF APPLIED PHYSICS 














It is advisable always to start by solving the 
difference equations with a very coarse net and then 
successively cut the net spacing in half. The 
previous solutions of the coarser nets are desir- 
able because at each halving of the net constant, 
the convergence is about four times as slow (cf. 
Fig. 2). Hence it is economical of time to get as 
good a function as possible from a coarse net to 
use as the initial trial function for the finer net, 
so as to keep the error in this function as small 
as possible. 

The net should be split up into as many nine- 
blocks as possible and the remainder into as many 
four-blocks as possible. Each nine-block or four- 








= 


Fic. 4. Location of 


~ 


interior points (circles) and boundary 
points (crosses). 


block the 
regular net spacing. These 12 or 8 surrounding 
points will ordinarily be interior points of the net 
but may, if a lattice line coincides with a section 


must be surrounded by points at 


of the boundary, be boundary points, as at the 
extreme right of Fig. 4. It seems advisable to 
observe these requirements rather than to go 
to the trouble of using formulas modified to take 
care of unusual distances of the points surround- 
ing the blocks.’ The more nine- and four-blocks 
one can use, the lower will be the eigenvalues and 
the faster the convergence. It occurs to one that 
it might be advantageous to use blocks which 
overlap, but we find that the extra labor involved 
with such blocks more than offsets the gain in 
rate of convergence. 

The reduction of the net spacing should be 
stopped when one has achieved the desired ac- 
curacy in the solution, which may of course not 
be greater than the accuracy with which the 
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boundary values are known. The amount by 
which the solution of the difference equations 
differs from that of the differential equation may 
be judged from the changes in final solutions for 
successively finer nets [cf. (2), which however 
may be expected to apply to the block procedures 
only on the average since the mode of division 
of successive nets into blocks will be different ]."! 


Choice of trial function 


The values of the function at the boundary points 
are of course determined by the given function 
which determines the boundary values for the 
region. In case this function is given only in the 
shape of a table of values at a discrete set of 
points around the boundary, the values at the 
boundary net-points must be determined by 
interpolation in this table, which is usually done 
most conveniently by plotting the given bound- 
ary values in some convenient way and reading 
off the values at the net points. 

The coarsest net used will usually converge so 
rapidly that one need not exercise particular care 
in guessing trial values for the interior points. 
When the net spacing is halved, the new trial 
function has the values of the final solution of 
the coarse net at the points which the two nets 
have in common and values interpolated from 
this solution at the remainder of the points. This 
interpolation may be very conveniently done by 
formulas which guarantee that Laplace's equation 
will be approximately satisfied at these points: 
If in Fig. 5 the black dots represent the points 
of the coarse net, the letters points at which 
interpolated values are necessary, one takes the 
value at a point marked a as the mean of the 
values at the four dots diagonally surrounding a, 
and then the value at a point marked 6 as the 


'! Some indication of the error in the difference solution 
may be obtained by comparing, at various points, the 
average of the four nearest neighbors and the average of 
the four diagonal neighbors. The error is certainly greater 
than the difference between these two averages. For the 
Liebmann procedure in a square net of pX p intervals, 
we can show that asymptotically the error is at most of 
the order of 0.06 p? times the difference between the 
value at a point and the average of its four diagonal 
neighbors. This is the case where the difference between 
the value at a point and the average of its four neighbors 
for the continuous solution is constant over the net. 
If the continuous solution is less regular the factor 0.06 is 
to be decreased. Such definite quantitative statements 
are not easy to make for other procedures. 
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mean of the values at the four nearest neighbors, 
two dots and two a points. 


a 


* b . 


Fic. 5. Interpolation. 


Improvement of points near boundary 


An interior point which has boundary values 
for neighbors will ordinarily be surrounded by 
two interior points at the regular net spacing h 


and two boundary points which are closer than 
h, say at distances sh and th. However it may be 
surrounded by three interior points, in which 
case either s or t may be taken as unity. The 
improved value of Yo in this case is chosen so that 
Laplace’s equation is approximately satisfied at 
the five points in question. The formula for this 
improved value, which is a special case of Eq. 
(23) of the Appendix, is 


st | 1 1 
+ : Vo 
l+s 

1 1 
+ —Wit Y; 
(1+) 1+¢ 


y=} 
(s+ols(1+s) 


(16) 


The complete value of the coefficients of y, and 
y, and the whole terms in y, and y, if these are 
constant boundary values should be calculated 
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once for all and written on the net. The value y, 
is then obtained very simply by multiplying y, 
and y, by known coefficients and adding a known 
constant. 


Axes of symmetry 


If there is an axis of symmetry such that the 
function values are known to be the same at 
corresponding points on opposite sides of this 
axis (AA might be such an axis in Fig. 4), it js 
only necessary to work with the net on one side 
of this axis. One imagines that each time a value 
is changed, the value at the corresponding point 
is simultaneously changed, and in improving 
values near or on the symmetry axis, one uses 
these imaginary values as needed. For example, 
if the left-hand row of points in Fig. 6 (p. 346) 
lie along a symmetry axis, and one is using the 
simple Liebmann formula for these points, the 
improvement formula for r is 


¥,=4(Wst2yit+y.). 


One obtains somewhat more rapid convergence 
in the Liebmann case if the net is chosen so that 
the symmetry axis is along one of the rows of net 
points rather than between two rows, and in the 
block procedures, if the blocks straddle the sym- 
metry axis rather than be located unsymmetri- 
cally. 

In a similar fashion one can handle a straight 
section of the boundary where the boundary 
condition is that the normal derivative of y have 
the value zero, or even have given nonzero 
values. 


Use of the difference function 


In most cases the labor of improvement will 
be materially reduced if one obtains the differ- 
ence between the first improved value and the 
original value, and then applies the improvement 
formula to this difference function. The difference 
function is easier to handle because it usually 
has many fewer significant figures than the 
original function, because it has zero boundary 
values, and because it permits of convenient 
extrapolation. In the notation of (5), (6), (7), 


y=ot+¢: 
the trial function equals the true solution plus 


JOURNAL OF APPLIED PHYSICS 





the error; when this is improved, one obtains 
Y=o4+T¢?. 


The difference between y and y we call the dif- 
ference function 


5=~—y=To-¢. 


Since this has zero boundary values, if we apply 
the improvement process to the difference func- 
tion we obtain successively the functions 


T5=T*d—T¢, 
T?5 = T'¢— T°4, 
T"-\§=T"¢—T"“\g. 


(17) 


If we add the sum of these functions to y¥, we 
obtain 


V+TS+T26+---+T"5=04T"G, (18) 


which is the same function we would have ob- 
tained by improving y itself times. 


Extrapolation 


The greatest advantage of using the difference 
function lies in the fact that it is possible at a 
certain point to make an accurate extrapolation 
to find the sum of the infinite series of terms of 
the type (17) which must be added to y to obtain 
the true solution w. [If m= « in (18), 


J+ETib=w, 
i=1 


since @ disappears after an infinite number of 
applications of T.] Since 6 has vanishing bound- 
ary values, it may be expanded in the series, 


cf. (13), 


6=)ip\ +b2p +. --+b,9™. (19) 


As T is applied successively to 6, the functions 
¢”, +++, @° disappear before ¢‘, and after a 
time T'6 will become a fairly smooth, pillow- 
shaped, one-signed function decreasing uniformly 
by a factor \; at each iteration. At this stage, one 
has 
T'6=co™, 
T*+16 =) ,co", 
T**?26=r2co™, 
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so that 
T*6+T* 6+ T7**26+ cee 
=cp(1+A,+A,2+---)=T*5 (1—A,), 


which furnishes the evaluation of the desired 
sum. The time saved by this extrapolation 
depends on the relative rate of convergence of 
the first and second eigenfunctions—after the 
second eigenfunction has been iterated out of 6, 
the first may be extrapolated out.” 


Resumé 


We conclude with a detailed resumé of the 
calculation for a given net, making use of the 
difference function and extrapolation. 

Having chosen the trial function, the improve- 
ment formula to be used for each point, and the 
order of going over the net, go over the trial 
function W once to obtain the function p. 

Calculate the difference function 6=~—y. This 
will have zero boundary values. It need not be 
computed at points whose unimproved function 
values are not used in the improvement of any 
other point. In the nine-block procedure, this 
usually applies to four of the nine points in each 
block, e.g. to points 41, 42, 50, 51 in the bottom 
center block of Fig. 1 if the blocks are improved 
in the order of reading a book. 

Improve the function 6 repeatedly in the same 
order and with the same improvement formulas 
as used for y, to obtain successively 76, T?6, 
T*5, ---. When this function has become fairly 
smooth, one-signed,'* and falling approximately 


2 For a large net in a boundary approximately square or 
circular, the extrapolation will reduce the number of 
necessary traverses by a factor of 2.5 since log \2/log A, = 2.5 
in this case, cf. (14c). This factor is reduced as the region 
becomes more rectangular or elongated, dropping to about 
1.3 for a rectangle 3 times as long as broad. But the factor 
is increased by the presence of axes of symmetry since in 
such cases the second eigenfunction will have a node along 
the axis of symmetry and the fact that all the successive 
approximation functions are automatically symmetrical 
means that the coefficient of this second eigenfunction in an 
expansion such as (19) is identically zero; in this case one 
will be ready for extrapolation as soon as the third eigen- 
function has disappeared. For a region approximately 
square or circular with two crossing symmetry axes, 
extrapolation reduces the labor by a factor of about 5. 
For a rectangular region with one axis of symmetry across 
the short way, the factor is a maximum of 3.2, falling to 
1.8 for a rectangle 3 times as long as broad. 

13 Cases sometimes occur where 7*6 is two-signed with an 
approximately stationary node and looks more like the 
second eigenfunction than the first. In this case it is 
profitable to extrapolate at this stage. 
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at a uniform rate, calculate the ratio T*6 T*—6 
at various points of the net. If the higher eigen- 
functions were completely eliminated, this ratio 
would be constant and equal to ,. It does not 
pay to wait for this to become strictly true before 
extrapolating. Fig. 2 and the 
footnote 6 assist in finding an early approxima- 
tion to the value of A. 

Using an approximate \, determined in this 
way, add together , Té, T*6, ---, T*%6, 
T*5 (1—X;) to obtain a new trial function y’. 
Apply the improvement process to y' to obtain yp’. 
(y’ need be determined only at those points at which 
its value is used in determining ’. In the case of 
the nine-block procedure this fact makes a con- 
siderable saving in the time of extrapolation.) The 
function y’ should be very close to the final 
function w, its departure from w is measured by 
the difference between y’ and y’. If necessary, 
calculate a new difference function 6'=y'—y' and 
repeat the procedure. The more constant one 
requires the ratio T*6 T*-'6 to be before extra- 
polating the better will be the function y’, but 
it seems advisable to make a fairly early extra- 
polation, and then continue the improvement, 


discussion of 


e . 


for two reasons. First, the new 6’ will in general 
be much smaller than T*é, and hence it will be 
much easier to improve it than to continue im- 
proving T*é. Second, there is a possibility that 
one may have made an error in obtaining p or 
one of the 7‘é’s, and if this has occurred, it will 
show up in y’, so it is advisable to make a fresh 
start before one does too much work, just as a 
check. 

In order to obtain a certain number of decimal 
places in the final solution w, it is necessary to 
carry more places in ¥ and 6 in order to avoid 
large rounding-off errors. It is necessary that a 
change of 1 (1—\A,) in the last figure carried in 6 


represent a change of less than 1 in the place 
where the last significant figure of w is to be 


located. It is desirable for ease of calculation of hy 


that the number of figures in 6 should be chosen 
more generously than this. For a net not exceed- 
ing a couple of hundred points, it is convenient to 
carry two more figures in W and 6 than one wants in 
the final solution. One then need carry the process 
only to the point where the estimated error of 
extrapolation is less than 50 in the last decimal 
place to get an error of less than } 
place of w. 


in the last 


With regard to the time and labor involved in 
the numerical solution of problems by these 
methods, no general statement can be made 
since so much depends upon the accuracy re- 
quired. To engineering accuracy of one percent 
a reasonably complicated problem such as arises 
in photoelasticity can be solved in two or three 
hours, after the boundary values have been 
determined.! For problems requiring greater 
accuracy we feel that the method is sufficiently 
convenient to be entirely practical and deserving 
of wide application to physical problems. 

We believe that methods similar to those dis- 
cussed above can be applied to the solution of 
Poisson’s equation Ay= p(x, y). We are planning 
to consider the numerical solution of this equa- 
tion in the near future." 

We are indebted to Professor L. H. Thomas for 
helpful discussions concerning approximation 
formulas. 


14 The heat conduction and vibration equation, Ay =)dy 
and Schrédinger’s equation Ay+p(x,y)¥=Ay are more 
complicated because an eigenvalue \ must be determined 
at the same time that y is determined. A generalization of 
the Liebmann procedure suitable for this type of equation 
has been considered by Kimball and Shortley, Phys. Rev. 
45, 815 (1934). 
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APPENDIX 


a. Approximation formulas 


We shall derive first the series which gives the value of 
a harmonic function ¥(x,y) at a given point (which we 
may conveniently take to be the origin 0,0) in terms of 
the function values yi, ¥2, +++, Yn at selected neighboring 
points Xi, Yi; X2, Ye; *; Xn, Yn. This will enable us to 
obtain improvement formulas for points of a net which 
are not surrounded by four regularly disposed points. In 
the neighborhood of the origin, any harmonic function can 
be represented by the series 


¥(x, vy) =aot+aur cos 06+ a2r? cos 20+--- 


+ ir sin 0+ 62r? sin 26+---, (20) 


where we have introduced polar coordinates r and 6; 
x=rcos 0, y=rsin@. The terms r” cos 6 and r” sin n@ 
here are homogeneous mth order polynomials in x and y, 
which are the real and imaginary parts of (x+7y)”. 

We now wish to approximate the value ¥(0,0)=yo of 
the arbitrary harmonic function (20) in terms of the 
function values ¥1, ¥2, -*:, Wn. We assume it to be the 
following linear combination 


n 


Yo= Db divi. 
al 


(21) 


When we substitute from (20) the values of these y’s, this 
equation becomes 
ay = a) Lb; +. Vir; cos 0; +2 Dir 2 cos 20;+--- 
+6; Dbir; sin 6:+B2Dbir;2 sin 20;+---. 
Since the a’s and #’s are to be arbitrary, we equate their 


coefficients to zero to obtain the following set of equations 
to determine the b’s of (21): 
2b; =1, 
Ddir; cos 0; = Db x; =(), 
Db,r; sin 6; = LDbiyi =(), 
Dbiri2 cos 26; = Db,(x2— yi?) =0, 


Lodi? sin 20; = Vbix:y; =0 


’ 


Of course, since the number of b’s is only n, one cannot 
satisfy this whole infinity of equations. One chooses the 
b’s to satisfy as many consecutive equations, beginning 
at the top, as possible; this corresponds to putting a 
harmonic function like (20) of as high an order as possible 
through the given points. 

It is Eq. (21) which is used as one of the set of difference 
equations which replace the Laplace differential equation, 
and which is also used for the improvement formula for 
the point in question. 

If n=5, we can in general satisfy the first five of equa- 
tions (22) and hence make use of a quadric interpolation 
function. The equation of footnote 9 is an example of 
this case. 

If n=4, we cannot in general find a set of b’s satisfying 
the first five equations of (22), so that we cannot use a 
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complete quadric interpolation formula. If however the 
four points are located on two straight lines passing through 
the origin and intersecting at 90°, there is a set of 6's which 
satisfy the first five equations of (22). We see this if we 
take these lines as axes of coordinates, in which case the 
fifth equation of (22) is automatically satisfied, and the 
same 6 values clearly satisfy the equations when each @, 
is decreased by a, corresponding to a rotation of coordinate 
axes through angle a. One of the coefficients a» or By» is 
undetermined by the function values at the four points, 
but all harmonic quadrics through the four points have 
the same value at the origin. For four points located at 
distances s,h, soh, ssh, ssh from the origin, as in the figure, 


(0, 5,h) 
s 
¥, 
(0,0) (34,0) 


*¥, °y, 


(-5,h, O) 


(0, - 
Pp. 3,4) 
%, 


the value of yo is given by 


( 1 4 1 ) 7 1 nN 1 
S1Se S354 adil 51(51 +52) vs So(Sy +S») v 


1 1 


+ Pere F S4(Sy+54) va 


(23) 
The improvement formula for points near a boundary (16) 
and points of a rectangular net with unequal spacing in the 
two directions follow immediately from this formula. For 
the special case s;=s2=5s;=5,4=1, where yo is just the 
average of the four neighboring points, the first 7 equa- 
tions of (22) are satisfied, so this formula is good to the 
third order—for all cubic harmonic polynomials, the 
value at a point is the mean of four equidistant neighbors. 

For n>5 we can in general satisfy Laplace’s equation 
to a higher order. For example if we take 8 points arranged 
in a square we find that 


x x x 
Ye Y2 Ys 


x * x 
Yo Yo Ya 


x” x x 


7 a Ye 


20f0=4 (Wit vetvsatys) tot voter tys. (24) 


The 6 values here satisfy the first fifteen equations of (22) 
so that Laplace’s equation is satisfied to the seventh order 
by this formula. 
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b. Change of net constant 


The considerations of the preceding section enable us 
to find suitable improvement formulas for a region such 
as that shown in Fig. 6 where the interval is being changed 


. . ° . . + 

° . . of . + 
. 7 . . . 2d oJ em + 

. . . . - eC et! 

eo . . 02 ed e/ ef 9 
we 2 « . . ef ek 

ou. . . . e e e 
. ° . . . . 

. . . . es ° Ph § 7 


Fic. 6. Change « 


net spacing. 


from the value / in one part of the net to double this 
value in another. Such a refinement of the net is some- 
times desirable near a region of the boundary where the 
function curves very rapidly. Typical formulas are 


¥a={Yat¥et¥;t+y¥.), 

¥n=4V t¥mtvotyo), 

Yo= 6 2votd stom t+ 2y,), 
Yie=b(Qvat3ity.+3yx), 

Vi= Pe (4yet+SyotdntyntS5y,). 


c. Accuracy of the difference solution. The seventh-order 
procedure 


We wish here to make a few general statements con- 
cerning the accuracy with which the converged solution 
obtained for a given net approximates to the continuous 
solution. In the Liebmann procedure the value of the 
function at a point is given in terms of the neighboring 
values correct to third-order terms. But the formula (24) 
satisfies the equation in the square surrounding each point 
to the seventh order. Hence in general one might expect 
that the use of (24) as an improvement formula will give 
a sufficiently closer approximation to the continuous 
solution than the other methods we have discussed to 
justify its use in spite of the extra labor entailed. The 
labor is greater because about three times the work per 
point is required in the improvement over that in the 
Liebmann or nine-block, while the rate of convergence is 
but slightly faster than that of the Liebmann and only 3 
as fast as the nine-block. (The asymptotic value for the 
seventh-order process, as we shall call this procedure, of 
— Pp? log \; is found to be ®5x* for a square net.) But in 
most engineering problems it is doubtful that one gains 
by using so accurate a formula, since it corresponds to 
using an interpolation formula which is more accurate 
than is indicated by the smoothness of the boundary 
values. Also, in the case of an irregular boundary, it is 
very inconvenient to attempt to use improvement formulas 
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which are more accurate than the second-order accuracy 
of (16) for points near the boundary. 

In cases where great accuracy is desired in the solution 
and is justified by the gift of a boundary function known 
analytically or to high accuracy and for which at least the 
fourth differences are smooth, the seventh-order formula 
would undoubtedly be very useful. To see how great a 
gain in accuracy is made by use of this formula in such a 
case, we chose an example in which the boundary values 
were taken so that the solution of the differential equation 
was logi 7* in the region with boundaries x = 100 and 110: 
y=-—5 and +5; net constant h=1 (81 interior points), 
In this case the maximum discrepancy from the con- 
tinuous solution was —514X 107!" for the Liebmann pro- 
cedure and only 700X10-' for the seventh-order pro- 
cedure. When the net constant was doubled the maximum 
error for the Liebmann procedure was — 1940 10-1, 
approximately four times the previous value, in accordance 
with (2). For the seventh-order procedure the error was 
46,000 10-5 with h=2, approximately 64 times the 
previous error, since with this more accurate interpolation 
formula the error varies as h® in place of h?. This is another 
great advantage of this procedure in cases where it can be 
used. 

On this same example, with h=1, the maximum error 
for the nine-block procedure was +108 X10~", and in 
each block there were both positive and negative errors 
while the Liebmann errors were everywhere larger and 
negative. These facts are explained by the fact that for 
log 7? in this region the value at a point is always greater 
than the average of its four neighbors, but less than the 
average of its four diagonal neighbors. The seventh-order 
formula (24) corrects for this by the proper mixture of 
averages taken in the two ways. The nine-block procedure 
apparently tends to make this correction on the average by 
calculating five of the points in the Liebmann way but the 
corner points diagonally. A similar type of correction of the 
error in the Liebmann procedure is possible in all cases, 
and it may be that in general the nine-block procedure will 
give a closer approximation than the Liebmann because of 
this partial correction. 


d. Proof of convergence of the iteration process 


Each of the single-point improvement formulas one uses 
is of the type (21) in which the improved value at a point 
is a weighted average (25; =1) of the values, some already 
improved, some unimproved, of the function at neighboring 
points. This is seen to imply that the improved value ¥; 
of (4) is in the last analysis a true weighted average of the 
values of y at all interior and boundary points, i.e. that 

ZTij+ DU ia=1. 
7 a 

The considerations of §2 furnish a proof of the con- 
vergence of the iteration process provided that we can 
prove that the magnitude of the greatest eigenvalue of T 
is less than unity. This will be proved if we can show that 
T applied to any function whatsoever of zero boundary 
values reduces the absolute maximum of the function. This 
is easily shown provided that all of the b's of the improve- 
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ment formula (21) are positive, since then all of the 7;,’s 
and Uva’s of (4) will be positive. Although it is conceiv- 
able that cases might arise in which negative values of b 
appear desirable, this is not the case in any situation which 
we have discussed or contemplated, and so seems unlikely. 
With positive b’s, the maximum value of a function is 
replaced upon improvement by a value which is a weighted 
mean (with all weights positive) of the values of the func- 
tion at all interior points and of the zero values of the 
function at the boundary points and this clearly reduces 
this maximum value.” No other point can get replaced by 
a value so large as this maximum, so the maximum of any 
function is definitely reduced upon operation with T. This 
is true in particular of the highest eigenfunction, so the 
highest eigenvalue must be less than unity. 


e. Calculation of asymptotic eigenvalues 


We shall illustrate the method of calculating the asymp- 
totic form of the eigenvalues as the number of net points 
increases by treating the case first discussed in the text 
in which the improvement is by blocks of v? points using 
such formulae that the difference equation (1) is satisfied 
at each point of the block by the improved values. [The 
eigenvalues for the more practical block procedures can be 
obtained in much the same way; we shall indicate in 
brackets the modifications of the calculation which are 
necessary in these cases. ] 

We shall consider first a square region containing m 
blocks on a side. Letter the points according to their 
location in the block, and number the blocks with two 
indices according to their position in the net, as in Fig. 7 
which applies to v=3. We denote the values of the eigen- 
function ¢ at the points aj;;, bj), -, simply by aj;, 
bij; +++; the values of the improved function ¢™ by 
a;;, bij, ---. If the blocks are improved in the order 11, 
71 3. -«« 3 82, 2a, Se, ; etc., the improved values in 
the ij block satisfy equations similar to the following, 
which are forv=3: 


—44;;+-fij+9:,+hij+l;=0, 
—4bij3+fijtljtea, j+&%, ;i=0, 
—4¢,+9:;+-fijt+di, ;-i.+b:, 
—4dij+hi,t+Giiteins, p+6%, 
—4@,, 41: ;+hij+bi, jnit4d, 

—4f,;+4;,+b);+0;+hi, 
—4g;,+4;,+¢:,4+di;4li, 
—4hi;+4;;4+dij,+6:,4+fi, 
—41;;44:;4+2:; +65 +9: 

[For other procedures one writes here the actual equations 


used in the improvement process, which are usually more 
complicated. ] 


The condition that @ represent an eigenfunction is that 


ai, =)ai;, bi, =)b;;, etc. 


(26) 


® For a function in which the maximum value is assumed at several 
neighboring points, it may be that for some order of improvement, 
not all of these maximum points will be reduced the first time over. 
They will all be reduced after several times over. The eigentunctions 


ire obviously not of this type, so this does not affect the convergence 
proof. 
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We shall not solve directly the problem of m? blocks 
with a row of zero values around the outside since these 
boundary conditions are difficult to express analytically. 
Rather we shall assume that the function value at a par- 
ticular type of point, say at the b points, is of the form 


ues . T . . Tv . - 
b;;= (1 +B :) sin ~ (t+x,) sin —(j+ »), (27) 
\ m? m m 


to terms of order 1/m?; here B, x,, and yp, are constants of 
the order of unity. Assuming a function of this form for 
each type of point amounts to replacing our problem by a 
periodic problem with “boundary" conditions a;, ;= 
;; etc. The 
choice (27) amounts to placing the “‘nodes”’ of the function 
in the blocks with 7 or 7=0, m, 2m, 
the function is of the order of 1 


value. For 


— Gi, j4m= —Gism, 73 Os, j= —Di, 54m = —Diam, 
+; in these blocks 
m times its maximum 
these homogeneous boundary conditions, 
(1—A,)m? will have the same asymptotic value as for the 


zero boundary conditions used in practice. 


Since A,;=1—O(1/m?) =1—k2?/m?+O(1/m?), (28) 


we need only satisfy the equations (25) to terms of the 
order of 1/m? in order to find the asymptotic value kx? of 
(1—Ai)m?. It turns out that it is possible to do this with 
functions of the type (27). We expand (27) in the form 


bite, i+e% = 114+LB— 1(x+ @)?— b (s+ We)? ]?/m?} 
sin ir/m sin jxr/m 
+ |(x,+ @)x/m} cosix/m sin jxr/m 
+ }(yo+ *&)r/m} sin ix/m cos jx/m 
+ | (xe + @ ) (n+ %&)2?/m?} cosin/m cos jr/m+O(1/m?), 


(29) 


with similar forms for dic@, jit, Cis@. ise: 

e =+1, 0, or —1; *=+1, 0, or —1. 
These expressions are to be substituted in the equations 

such as (25), and the coefficients of the sin sin, sin cos, 


etc. Here 


cos sin, and cos cos terms to be equated to zero. From the 
form of (25) and the coefficients of cos sin in (29) we see 
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that the x+ @ values must increase in uniform steps as 
we move to the right, so that, except for an additive 
constant, 


tq= 1 


r= 0, 
as Xg _ Xe = 5 ’ 


=x tl=m+1=+3 


xe+1 


for v=3. In general the corresponding values will be 
(—v—1)/2v, (—v+1)/2v, +--+, (v—1)/20, (v+1)/22, 
arranged similarly in the blocks. The y values are the same 
as the x values with the block turned through 90°. With 
these values it is easily verified that the cos cos equations 
are automatically satisfied. [For all other procedures for 
which we have calculated asymptotic values, these same 
x and y values have satisfied the cos sin, sin cos, and cos cos 
equations. ] 

Finally, equating the coefficients of the sin sin terms 
furnishes a set of v linear equations to determine the v* 
constants A, B, C, -, and also k. We see that each of 
the capital letters occurs —4 times in one equation and 
once in four equations, so that if the whole v? equations 
be added A, B, C, --- drop out leaving a single equation 
for k. [The equations obtained from setting the coefficients 
of sin sin equal to zero are in general unchanged by adding 
the same constant to A, B, C, 
multiplication of ¢ by a constant of the magnitude of 
1+0(1/m?) 
between them. In the twenty-five block procedure it was 


—this corresponds to 
so we can assume one arbitrary relation 


to assume the sum of the values at the four a 
points of Fig. 3 to be zero before a linear combination of 


necessary 


the equations could be found which eliminated all unknowns 
but &. In the other cases we have considered, a satisfactory 
linear combination could be made without such an assump- 
tion. ] From (26) we see that this 
equation once for each barred value so k& occurs once for 


—kr*/m*? occurs in 


each unbarred value. Since unbarred values occur only 
on two edges of the block, the coefficient of kx?/m? in the 


resultant equation is 2v. Finally we must calculate the 
constant term. Since each letter occurs —4 times and +4 
times in (25), the x? and y? parts of (29) cancel, leaving 
only the terms involving @ and *. The @ terms involve 
+x for each of the v points bordering the ij block to the 
left and —x for the points bordering to the right; giving 
in particular for 0=3, xa+Xp+x-—X.—X1—%}=6-}=2; 
in general, 2v-(v—1)/2v=(v—1) in units 2*/m?. One has 
also —}@? for each of these points, or —v altogether, 
The whole contribution of the @ terms is then —1; the 
% terms also give —1, total —2z2?/m?. 

The equation resulting upon addition of all the equations 
of the type (25) is then 


_ 2 / yy? — 22/2? = 9772 
or 1—)A, =kr?/m? = x? /vm? = v2" / p*, 


where p=vm (asymptotically) is the number of intervals 
on a side. This is the result discussed in the text. 

For a rectangular region one obtains a value equal to 
(14b) with vz? written for 7?. One can derive this by the 
method just outlined if one takes m blocks in one direction, 
l blocks in the other. Asymptotic values for the other 
eigenfunctions can be shown to be given by (14c) with vz? 
written for x? by taking in place of (27) a product of 
sin (ux/m)(i+x,) sin (vr/l)(j+y). [For any other pro- 
cedure, (14a), (14b) and (14c) are also obtained, with 
const. Xz? replacing 7. ] 

Finally we might remark that if in the improvement of 
a block the original @ values are used for all boundary 
values, in place of the improved value @ where they are 
available ; e.g. if in (25) the bars are removed from letters 
with subscripts i—1, j and 7, j—1; the asymptotic values 
of 1—X come out only half as great as before. This is true 
of any procedure in which the points of a block are treated 
symmetrically. Hence we see that the device of using 
improved values as fast as they are available in the im- 
provement of other points just doubles the rate of con- 
vergence over a similar procedure in which improved values 
are not used. 
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